AFRL-VS-HA-TR-2006-1047 


Simultaneous  Inversion  of  Receiver  Functions,  Multi- 
Mode  Dispersion,  and  Travel-Time  Tomography  for 
Lithospheric  Structure  Beneath  the  Middle  East  and 
North  Africa 


Charles  J.  Ammon 
Minoo  Kosarian 
Robert  B.  Hermann 


Saint  Louis  University 
Office  of  Research  Services 
3634  Lindeil  Blvd. 

Saint  Louis,  MO  63108-3342 


Final  Report 


1  February  2006 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


AIR  FORCE  RESEARCH  LABOR  ATORY 


Space  Vehicles  Directorate 
29  Randolph  Road 


\  AIR  FORCE  MATERIEL  COMMAND 
/  Hanscom  AFB,  MA  01731-3010 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
AFRL- VS-HA-TR-2006- 1 047 


/signed/ 


/signed/ 


ROBERT  RAISTRICK  ROBERT  BELAND,  Chief 

Contract  Manager  Battlespace  Surveillance  Innovation  Center 


This  report  has  been  reviewed  by  the  ESC  Public  Affairs  Office  (PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS). 

Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical 
Information  Center  (DTIC).  All  others  should  apply  to  the  National  Technical 
Information  Service. 

If  your  address  has  changed,  if  you  wish  to  be  removed  from  the  mailing  list,  or  if 
the  addressee  is  no  longer  employed  by  your  organization,  please  notify 
AFRL/VSEM,  29  Randolph  Rd.,  Hanscom  AFB,  MA  01731-3010.  This  will  assist 
us  in  maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a 
specific  document  require  that  it  be  returned. 

Using  Government  drawings,  specifications,  or  other  data  included  in  this  document 
for  any  purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the 
U  S.  Government.  The  fact  that  the  Government  formulated  or  supplied  the  drawings, 
specifications,  or  other  data  does  not  license  the  holder  or  any  other  person  or 
corporation;  or  convey  any  rights  or  permission  to  manufacture,  use,  or  sell  any 
patented  invention  that  may  relate  to  them. 

This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange 
and  its  publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its 
ideas  or  findings. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 

data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202' 

4302  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

01-02-2006  |  Final  Report 

3.  DATES  COVERED  (From  -  To) 

1  May  2002  -  31  Jan  2006 

4  TITLE  AND  SUBTITLE 

Simultaneous  Inversion  of  Receiver  Functions,  Multi-Mode  Dispersion,  and  Travel-Time 
Tomography  for  Lithospheric  Structure  Beneath  the  Middle  East  and  North  Africa 

5a.  CONTRACT  NUMBER 

DTRA0 1 -02-C-0038 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

S49012 

6.  AUTHOR(S) 

Charles  J.  Ammon*,  Minoo  Kosarian*  and  Robert  B.  Hermann  * 

5d.  PROJECT  NUMBER 

DTRA 

5e.  TASK  NUMBER 

OT 

5f.  WORK  UNIT  NUMBER 

A1 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

St.  Louis  University 

Office  of  Research  Services 

3634  Lindell  Blvd. 

Saint  Louis,  MO  63108-3342 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFRL/VSBYE 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

AFRL-VS-HA-TR-2006-1047  | 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  Unlimited. 


1  Penn  State  University 

2  Saint  Louis  University 


3  Lawrence  Livermore  National  Laboratory 

4  Los  Alamos  National  Laboratory 


13.  SUPPLEMENTARY  NOTES 

*  Penn  State  University 
**  Saint  Louis  University 


14.  ABSTRACT 

An  important  component  to  understanding  the  evolution  of  the  continental  lithosphere  is  to  improve  our  knowledge  on  lower  continental  composition. 
Contribution  towards  this  goal,  we  perform  receiver  function  analysis  using  teleseismic  waveforms  recorded  at  permanent  and  temporary  broadband 
seismic  station  slocated  in  Middle  East,  Europe,  Asia,  and  North  Africa.  Two  hundred  and  twenty  six  stations  recording  a  total  of  about  6,000  teleseicmic 
events  producing  more  than  100,000  seismograms  have  been  investigated.  The  distribution  includes  72  stations  in  the  Middle  East,  57  stations  in  Europe, 
60  stations  in  Asia,  and  37  stations  in  central  and  north  Africa.  We  have  examined  receiver  functions  for  213  of  stations  (best  data)  in  the  period  of  1990  - 
2004  and  applied  the  receiver  function  stacking  procedure  of  Zhuand  Kanamori  [2000]  to  estimate  Poisson’s  ratio  and  crustal  thickness.  We  have  divided 
the  research  area  according  to  five  tectonics  environments,  explicitly  Shields,  Platform,  Paleozoic  orogenic  belts,  MesozoicCenozoic  orogenic  belts,  and 
rift  zones  based  on  Condie’s  [1989]  simplified  classifications.  The  results  from  this  study  shows  lower  value  of  Poisson’s  ratio  o  =  0.25  for  Shield  and 
Platform  compare  to  the  Orogenic- belts  with  a  =  0.27.  Crustal  thickness  for  Shield  and  Platform  show  the  value  of  38  km  and  43  km  respectively,  while 
for  the  Orogenic  belts  we  found  a  value  of  37  km  for  Paleozoic  belts  and  39  km  for  Mesozoic-Cenozoic  belts,  although  the  range  of  thicknesses  for  the 
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Chapter  1 

Introduction  &c  Method  Overview 


Of  the  many  successes  in  recent  seismic  verification  research,  significant  advances 
in  the  application  of  kriging  to  construct  travel  time  and  seismic  phase  amplitude 
correction  surfaces  from  available  sparse  calibration  points  (Schultz  et  al,  1998; 
Phillips,  1999;  Phillips  et  al  1998;  Rodgers  et  al  1999;  Pasyanos,  2000)  are  perhaps 
most  notable.  Information  gleaned  from  well  placed  calibration  events  significantly 
reduces  mislocation  to  the  extent  that  in  some  instances  may  remove  the  need  for 
regional  crustal  models  to  improve  seismic  event  locations.  Of  course  the  success 
of  location  and  event-size  methods  relying  on  kriged  correction  surfaces  requires 
calibration  events  with  ground-truth  location  information.  Not  all  places  of  inter¬ 
est  have  sufficient  events  or  ground  truth  information  and  some  extrapolation  is 
required  to  incorporate  the  knowledge  into  standard  processing  and  analysis  proce¬ 
dures.  In  these  regions,  basic  investigations  into  the  character  of  crustal  propagation 
are  required  for  several  reasons.  First  to  help  improve  event  location,  and  second 
to  prepare  for  and  aid  the  ultimate  goal  of  event  identification.  If  a  unique  event  is 
observed,  significant  resources  are  directed  toward  a  timely  identification.  Implica¬ 
tions  for  national  security  dictate  that  the  identification  be  robust,  which  requires 
that  we  exploit  all  available  observations  and  tools. 

One  of  the  principal  problems  in  identifying  and  locating  small  seismic  events  is 
the  low-amplitude  of  observed  signals.  Without  data  from  a  nearby  station,  we  are 
forced  to  deal  with  seismic  observations  that  are  small  in  amplitude  and  composed 
of  short-period  signals  that  are  strongly  distorted  by  propagation.  At  a  few  thou¬ 
sand  kilometers  from  the  source,  the  largest  signals  generated  by  a  shallow  event 
are  short-period  surface  waves.  To  reliably  identify  and  measure  these  signals  some 
waveform  processing  such  as  phase-match  filtering  (Herrin  and  Goforth,  1977)  will 
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Figure  1.1:  Global  group  velocity  variations  from  the  Harvard  global  mapping 
project  [Larson  and  Ekstrom,  2001].  Although  the  period  range  is  too  long  for 
the  study  of  small  events  (which  have  little  useful  long-period  energy),  the  models 
provide  valuable  information  on  structure  for  regional  studies  of  the  lithosphere  that 
increase  resolution  using  shorter  period  observations. 


be  necessary.  The  construction  of  optimal  phase-match  filters  requires  knowledge  of 
the  dispersive  character  of  the  upper  crust  between  the  source  and  the  seismometer 
and  that  knowledge  can  be  derived  from  either  direct  measurements  of  the  dispersion 
of  calibration  signals  (when  available)  or  indirect  knowledge  of  the  shear  velocity  of 
the  shallow  Earth  along  the  wave  path.  Undoubtedly  the  best  approach  to  estimat¬ 
ing  surface- wave  dispersion  along  a  specific  path  is  direct  observation.  Calibration 
events  with  ground-truth  location  information  can  provide  region-specific  dispersion 
characteristics  for  areas  of  interest  that  have  well-studied  seismic  activity.  Unfortu¬ 
nately  not  all  places  of  interest  have  frequent  events  or  ground  truth  information  and 
some  extrapolation  is  required  to  incorporate  the  knowledge  in  standard  processing 
and  analysis  procedures. 

Confident  and  reliable  characterization  of  regional  structure  requires  effort.  With 
sufficient  observations  of  paths  crossing  the  area  of  interest  we  can  use  tomographic 
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Figure  1.2:  Illustration  of  the  variation  of  signal  frequency  content  for  a  source- 
receiver  distance  of  about  1000  km.  The  signals  on  the  left  were  generated  by 
earthquakes  in  the  Tien  Shan  region  and  recorded  at  station  WMQ  in  central  Asia. 
Those  on  the  right  are  band-pass  filtered  version  of  the  raw  signals  with  a  filter  pass¬ 
ing  periods  between  25  and  15  second  (spanning  the  20s  Ms  measurement  range). 
The  surface  wave  of  the  Mw  3.6  event  is  visible  in  the  noisy  trace  at  the  bottom.  En- 
hacning  the  signal  to  make  an  Ms  measurement  would  require  an  appropriate  short- 
period  phase-match  filter.  Mw  values  are  from  moment-tensor  inversion  (Ghose  et 
ah,  1998). 
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imaging  methods  to  construct  a  smooth  version  of  the  regional  dispersion  structure. 
Much  effort,  and  much  success  has  been  achieved  in  efforts  to  map  surface-wave 
dispersion  variations  across  the  planet  (Figure  1.1).  This  approach  has  shown  great 
promise  as  researchers  apply  the  method  to  intermediate  and  short-period  surface- 
wave  observations,  but  the  task  becomes  more  difficult  as  one  proceeds  to  periods 
less  than  15  to  20  seconds.  These  periods,  which  are  critical  to  enhancing  signals 
from  small  events  (Figure  2),  are  difficult  to  observe  routinely  and  Rayleigh- wave 
dispersion  tomograms  for  the  short  periods  (and  the  longest  periods)  often  have  the 
least  number  of  observations  in  tomographic  analyses.  Compounding  the  problem 
is  the  fact  that  short-period  surface  waves  are  sensitive  to  structures  that  are  small 
enough  to  require  more  paths  to  image  accurately.  One  approach  to  constructing 
the  necessary  knowledge  is  to  combine  the  longer-period  tomography  observations 
with  other  seismic  data  to  improve  regional  shear  velocity  maps  which  can  then  be 
used  to  predict  path-dependent,  short-period,  surface-wave  dispersion.  That  is  the 
avenue  we  propose  to  pursue  with  the  goal  of  constructing  site-specific  structure 
map  of  the  crust  and  upper  mantle  beneath  China. 


Composite  Seismic  Data  Inversions 


Seismic  studies  of  earth  structure  have  historically  proceeded  with  a  divide  and  con¬ 
quer  approach.  Different  parts  of  seismograms  were  isolated  and  modeled  by  differ¬ 
ent  methodologies  and  often  different  researchers.  The  divide  and  specialize  strategy 
has  computational  advantages  that  are  much  less  important  today,  but  our  experi¬ 
ence  in  combining  historically  disjoint  observations  remains  limited.  Subsurface  ge¬ 
ology  generally  has  a  broad  wavenumber  spectrum  (Figure  2),  containing  sharp,  or 
high- wavenumber,  changes  in  velocity  near  Earths  surface,  at  the  sediment-basement 
transition,  near  the  crust-mantle  boundary  (usually),  and  in  the  upper-mantle  tran¬ 
sition  zone,  and  smooth  low-wavenumber  variations  in  regions  of  relatively  uniform 
geologic  structure.  Access  to  the  full  spectrum  of  earth  structure  requires  that  we 
exploit  signals  that  span  a  wide  frequency  range  and  that  are  sensitive  to  the  entire 
spectrum  of  heterogeneity.  For  example,  surface-waves,  travel  times,  and  direct- 
wave  amplitudes  are  sensitive  to  smooth  variations  in  earth  structure;  reflected  and 
converted  waves  are  sensitive  to  the  velocity  contrasts.  Combining  seismic  data 
in  joint  inversions  is  an  obvious  approach  to  improve  estimates  of  earth  structure. 
Refraction  seismologists  have  long  used  seismic  wave  travel  times,  amplitudes,  and 
gravity  variations  to  study  the  structure  along  one-dimensional  surface  profiles  that 
illuminate  the  two-dimensional  subsurface  geology.  To  successfully  combine  data  in 
an  inversion,  we  must  insure  that  all  the  data  are  sensitive  to  the  same  (or  related) 
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Figure  1.3:  Each  triangle  represents  a  seismic  station  for  which  we  analyzed  receiver 
functions.  Not  all  stations  had  structures  suitable  for  a  joint  inversion  of  receiver 
functions  and  surface-wave  dispersion;  some  receiver  functions  indicated  strong  het¬ 
erogeneity,  for  which  simple  plane-layered  modeling  is  inappropriate. 
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physical  quantities  and  that  they  sample  or  average  structure  over  comparable  length 
scales.  Recent  advances  in  surface-wave  tomography  have  provided  an  opportunity 
to  combine  localized  surface-wave  dispersion  estimates  with  other  data  such  as  P- 
and  S-wave  receiver  functions. 

No  one  doubts  that  using  all  available  data  to  construct  earth  models  is  the  ap¬ 
propriate  direction  for  future  studies,  but  few  careful  studies  of  the  problems  and 
specific  benefits  of  joint  inversions  are  documented  in  the  literature.  Some  recent 
studies  have  included  modeling  complete  waveform  (usually  intermediate-to-long 
period  body  and  surface  waves)  to  extract  source  or  path  structure  information. 
Although  waveform  inversion  is  simple  and  straight  forward,  it  is  unclear  if  it  is 
the  best  way  to  extract  information  about  earth  structure.  Unlike  a  surface-wave 
dispersion  study  where  the  observations  at  a  particular  period  are  equally  weighted, 
the  period-dependent  weighting  in  a  waveform  is  not  even  and  simple  norms  such  as 
LI  and  L2  misfit  measures  are  easily  dominated  by  Airy  phases.  Such  observations 
are  also  not  completely  independent  since  both  the  body  and  surface-wave  group 
delays  can  be  affected  by  mislocation  and  origin  time  errors.  Finally,  we  note  one 
other  reason  for  utilizing  dispersion  measurements.  The  mechanism  by  which  dif¬ 
ferent  seismic  waves  traveling  through  Earth  average  heterogeneity  is  an  area  that 
requires  some  investigation.  Does  a  one-dimensional  parameterization  of  subsurface 
structure  contain  biases  resulting  from  an  uneven  averaging  of  the  structure  by  the 
seismic  wave  field?  To  avoid  such  issues  we  chose  to  invert  tomographically  sepa¬ 
rated  variations  in  dispersion  parameters  and  receiver  functions,  which  sample  the 
structure  within  a  few  10s  of  kilometers  of  the  seismic  recording  station.  This  in¬ 
sures  a  relatively  local  sampling  to  the  station  and  in  some  instances  better  justifies 
the  simple  one-dimensional  parameterization. 

VVe  have  successfully  applied  the  joint  inversion  of  receiver  functions  on  data  from 
across  western  Eurasia  and  central  and  northern  Africa  (Figure  1.3).  We  illus¬ 
trate  the  recent  results  of  our  efforts  at  joint  receiver-function  fundamental  surface- 
wave  inversion  later,  but  our  work  has  produced  two  important  results  that  de¬ 
serve  mention  before  we  elaborate  into  a  detailed  discussion.  First,  the  addition 
of  intermediate-period  group  velocities  (10  to  50  seconds  period  band)  to  receiver 
functions  greatly  reduces  trade  off  between  depth  and  velocity  above  velocity  con¬ 
trasts  in  receiver-function  models  by  adding  constraints  on  broad  averages  of  the 
crustal  shear  velocity  (receiver  functions  provide  solid  information  on  the  thickness 
and  magnitude  of  the  major  velocity  transitions  in  the  crust).  Second,  initial  appli¬ 
cations  (in  regions  of  relatively  simple  structure)  have  produced  models  that  agree 
well  with  existing  estimates  of  structure  from  refraction  profiling. 
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Figure  1.4:  Expansion  of  a  shear- velocity  model  (dark  line  in  top  panel)  in  terms  of 
Chebyshev  Coefficients.  The  top  panel  shows  the  sum  of  two,  five  and  fifty  terms 
respectively.  The  lower  panels  show  the  model  spectra,  that  on  the  right  focuses  on 
the  smaller  amplitude,  higher  wavenumber  coefficients.  In  practice,  the  two  data  sets 
are  band  limited  and  do  not  cover  the  entire  spectrum  of  the  model  -  the  illustration 
is  conceptual. 
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Receiver  Functions 


Receiver  functions  represent  the  response  of  a  location  structure  to  an  incoming 
plane  wave  (e.g.  Phinney,  1964;  Langston,  1979).  The  response  consists  of  the  di¬ 
rect  P-wave  arrival  followed  by  smaller  amplitude  P  and  S  waves  which  reverberate 
in  the  local  structure.  By  modeling  the  arrival  time  and  amplitude  of  the  rever¬ 
berations  we  can  extract  information  on  the  structure  surrounding  a  seismometer 
with  a  radius  of  about  40  km  (Owens  et  al.  1984).  The  use  of  teleseismic  P  waves 
makes  the  technique  suitable  for  temporary  deployments  since  sufficient  observa¬ 
tions  accrue  quickly  (provided  earthquake  sources  in  the  distance  range  of  30  to  90 
degrees  are  available,  which  is  true  of  most  locations).  To  isolate  the  near-receiver 
response,  the  vertical  component  of  ground  motion  is  deconvolved  from  the  hori¬ 
zontals,  removing  the  signature  of  the  source  and  near-source  structure  (Langston, 
1979).  The  cost  of  the  deconvolution  procedure  is  the  loss  of  P-wave  multiples  in 
the  resulting  time  series  which  is  easily  accounted  for  in  synthetic  waveforms  and 
results  in  an  enhanced  sensitivity  to  shear- velocity  (Langston,  1979;  Owens  et  al., 
1984;  Ammon,  1991).  Receiver  functions  are  primarily  sensitive  to  the  sharpness 
of  velocity  contrasts  beneath  the  receiver  and  the  vertical  travel  time  between  the 
contrasts  (Ammon  et  al.,  1990),  and  receiver  function  waveforms  are  relatively  easy 
to  model.  Owens  et  al  (1984)  formulated  the  first  time-domain  receiver  function 
inversion  and  Ammon  et  al.  (1990)  extended  the  procedure  to  include  reflectiv¬ 
ity  waveform  modeling,  smoothness  constraints,  and  investigations  of  initial  model 
dependence  on  the  results. 


Receiver  Function  Estimation.  When  the  data  are  high-quality  and  the  re¬ 
ceiver  structure  is  not  too  complex,  the  choice  of  a  deconvolution  procedure  does  not 
make  much  difference.  However,  when  the  noise  in  the  seismograms  is  substantial,  or 
the  receiver  structure  is  complex,  different  deconvolution  approaches  have  strengths 
and  weaknesses.  We  computed  receiver  functions  using  the  iterative  time-domain 
deconvolution  procedure  described  by  Ligorra  and  Ammon  (1999).  We  prefer  the 
iterative  approach,  which  is  based  on  the  Kikuchi  and  Kanamori  (1982)  source-time 
function  estimation  algorithm,  for  several  reasons.  First,  in  the  iterative  approach 
the  receiver  function  is  constructed  by  a  sum  of  Gaussian  pulses  which  produces  a 
flat  spectrum  at  the  longest  periods.  The  flat  long-period  spectrum  can  be  viewed 
as  a  priori  information  that  helps  reduce  side-lobes  that  may  result  of  spectral  or 
singular-value  truncation  stabilization  procedures.  The  reduction  of  side-lobes  eases 
the  interpretation  and  helps  stabilize  low-frequency  receiver  functions.  Second,  the 
iterative  approach  constructs  a  causal  receiver  function,  which  is  what  we  expect 
in  all  cases  of  reasonable  earth  structure.  This  is  a  subtle  difference  from  spectral 
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techniques  (e.g.  Langston,  1979;  Park  and  Levin,  2000)  which  can  always  introduce 
a  component  to  the  signal  before  the  P-wave.  The  acausal  component  of  the  spec¬ 
tral  signal  may  be  small,  but  still  important  to  the  satisfaction  of  the  convolutional 
model  that  defines  a  receiver  function 

(IT)  R(t)  =  Z(t)*ER(t). 

where  R(t)  and  Z(t)  are  the  radial  and  vertical  components  of  ground  motion,  and 
Efi(t)  is  the  radial  receiver  function  (a  similar  equation  holds  for  the  transverse 
component).  The  point  is  that  even  when  the  receiver  function  estimation  is  unsta¬ 
ble,  spectral  deconvolutions  may  satisfy  (1.1)  quite  well.  The  iterative  time-domain 
approach,  which  can  be  restricted  to  produce  the  best  causal  solution,  may  not  al¬ 
ways  satisfy  (1.1).  Experienced  modelers  have  always  been  able  to  identify  failed 
receiver  functions,  but  the  misfit  to  (1.1)  available  from  iterative  deconvolutions 
provides  quantitative  information  that  can  be  used  when  stacking  signals,  or  in  ex¬ 
treme  cases,  to  discard  obviously  failed  deconvolutions.  In  our  case  we  find  using 
a  threshold  cut-off  of  80-90%  of  the  radial  power  fit  allows  us  to  quickly  discard 
poorly  constrained  deconvolution  results,  enabling  an  efficient  and  objective  select 
of  the  data  to  include  in  further  analysis.  We  generally  allow  the  iterative  inversion 
to  include  more  than  600  arrivals  in  the  receiver-function  estimate.  For  longer  time 
windows,  we  allow  up  to  1200.  The  important  point  is  not  to  under  sample  the 
response  by  including  a  limited  number  of  iterations  in  the  deconvolution  (a  sign  of 
under-sampling  are  flat  regions  in  the  receiver- function  estimate). 

The  receiver  function  data  flow  used  in  our  analysis  is  summarized  in  Figure  1.5. 
The  original  P- waveforms  are  obtained  from  the  IRIS  DMC  for  all  events  with  mag¬ 
nitudes  greater  than  5.5.  Each  three  component  P-waveform  is  visually  examined 
and  the  P-wave  arrivals  are  marked.  Events  with  good  visual  signal-to-noise  ratios 
are  used  with  the  iterative  deconvolution  method  outlined  in  Ligorria  and  Ammon 
[1999].  Those  that  have  convolution  misfits  less  than  20%  are  included  in  later 
analyses  that  include  receiver-function  stacking  and  shear-  velocity  inversion. 


Poissons  Ratio  &;  Crustal  Thickness  Estimation.  As  a  first  step  in  the  re¬ 
ceiver  function  analysis  we  use  the  receiver  function  stacking  method  of  Zhu  and 
Kanamori  (2000)  to  estimate  the  crustal  thickness  and  Vp/Vs  velocity  ratio  (or 
Poissons  ratio).  The  stacking  method  makes  a  rather  limiting  assumption  of  a  uni¬ 
form  crust  but  the  analysis  provides  good  estimates  of  these  quantities  when  the 
structure  is  relatively  simple.  The  estimated  values  of  Poissons  ratios  can  be  used 
in  subsequent  inversions  which  require  some  assumed  value  of  bulk  crustal  Poissons 
ratio.  We  summarize  our  results  in  Figure  3.  Although  it  is  difficult  to  just  cluster 
regions  of  the  continent  into  simple  classifications,  we  follow  Zandt  and  Ammon 
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Figure  1.5:  Sketch  of  the  data  processing  flow  for  the  receiver  function  estimation 
used  in  the  project. 


(1995)  and  separate  our  data  into  tectonic-age  based  groups  with  the  exception  of 
Tibet,  which  we  remove  keep  as  a  separate  subset  of  the  Mesozoic-Cenozoic  regions. 
Our  results  are  similar  to  those  of  Zandt  and  Ammon  (1995)  with  the  exception 
that  our  larger  survey  of  shields  (including  Archean  shields)  does  not  show  an  un¬ 
usually  high  Poissons  ratio  as  was  suggested  in  the  earlier  data  set.  Only  7  stations 
are  common  to  both  studies,  so  the  differences  are  a  sampling  result.  Interestingly, 
we  see  no  differences  between  Archean  and  Proterozoic  shields.  Crustal  thickness 
estimates  agree  on  average  with  the  global  crustal  model  2.0,  but  at  times  the  differ¬ 
ences  are  significant  (greater  than  5  km).  The  mean  difference  between  our  results 
and  Crust  2.0  is  about  1  km  and  the  standard  deviation  of  the  differences  is  about 
7  km,  which  is  larger  than  the  resolution  we  have  on  stations  with  stable  receiver 
functions  (probably  between,  2.5  to  5  km).  The  numbers  agree  more  closely  when 
we  rank  our  estimates  using  the  complexity  of  the  observed  receiver  functions.  At 
the  bottom  of  Figure  3  we  show  the  estimated  Vp/Vs  values  for  each  station.  The 
shaded  region  shows  a  liberal  range  that  is  typical  for  the  continental  crust.  For  the 
most  part  our  results  are  showing  up  lower  than  average  (Zandt  and  Ammon,  1995). 


Fundamental-Mode  Surface- Wave  Dispersion 


In  the  time  domain,  surface  waves  often  have  uneven  period-dependent  weighting 
due  to  the  existence  of  Airy  phases  and  differences  in  source  excitation.  Dispersion 
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measurements  equalize  the  weighting  and  have  been  used  to  estimate  earth  structure 
for  more  than  four  decades  (e.g.  Haskell,  1954,  Brune,  1969;  Der,  et  ah,  1970; 
Braile  and  Keller.  1975).  Observed  dispersion  can  be  used  to  constrain  portions  of 
the  subsurface  structure  along  the  propagation  path  but  substantial  trade-off  in  the 
details  of  the  lower  crust  and  upper  mantle  are  difficult  to  control  without  additional 
information.  Further,  dispersion  measurements  represent  an  average  between  the 
source  and  the  receiver,  and  thus  average  the  structure  over  a  potentially  large 
path.  To  estimate  local  variation  near  a  receiver  we  will  rely  on  dispersion  estimates 
from  tomographic  analyses  that  geographically  localize  dispersion  observations.  The 
mathematics  of  the  problem  are  quite  simple.  In  the  same  way  a  linearized  receiver 
function  is  formulated  to  iteratively  estimate  the  variations  in  velocity  beneath  the 
receiver,  the  dispersion  measurements  are  used  to  iteratively  solve  the  nonlinear 
problem  relating  the  observed  frequency-dependent  velocities  on  to  the  shear-wave 
velocity  along  the  path.  Partial  derivatives  are  readily  computed  using  perturbation 
theory  (e.g.  Rodi  et  al.,  1975). 

A  companion  collaboration  of  this  work  was  the  continued  improvement  of  the  group- 
velocity  dispersion  models  for  the  MENA  region  by  researchers  at  LLNL  [Pasyanos 
and  Walter,  1999].  We  updated  the  dispersion  models  three  times  during  our  anal¬ 
yses,  as  the  information  was  provided  by  LLNL.  We  focus  little  discussion  in  this 
report  on  the  model  since  the  work  was  done  elsewhere  by  others.  However,  our 
work  would  not  have  been  possible  without  the  related  effort  at  LLNL  and  earlier 
studies  by  groups  at  Harvard  [Larson  and  Ekstrom,  2000]. 


Shear- Velocity  Inversion  with  Constraints 


The  receiver  function  is  sensitive  to  velocity  transitions  and  vertical  travel  times, 
surface-wave  dispersion  measurements  are  sensitive  to  aver-ages  of  the  velocities, 
and  relatively  insensitive  to  sharp  velocity  contrasts.  The  complementary  nature 
of  the  signals  makes  them  ideal  selections  for  joint  study  because  they  can  fill  in 
resolution  gaps  of  each  data  set.  Ammon  and  Zandt  (1993)  pointed  this  out  in  a 
study  of  the  Landers  region  of  southern  California  (although  for  their  specific  case, 
available  observations  were  unsuitable  to  resolve  subtle  features  in  the  lower  crust) 
and  Ozalaybey  et  al.  (1997)  and  Last  et  al.  (1997)  have  performed  complementary 
analyses  of  surface-wave  dispersion  and  receiver  functions.  The  mechanics  of  the 
inversion  are  relatively  simple  since  partial  derivatives  of  dispersion  observations 
(Herrmann,  1995)  and  receiver  functions  waveforms  (e.g.,  Randall.  1989,  Ammon 
et.  al,  1990)  can  be  calculated  quickly  and  accurately. 
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The  inversion  is  formulated  to  allow  convenient  adjustment  of  the  relative  impor¬ 
tance  of  the  different  signals  on  the  resulting  models.  We  use  an  iterative  least- 
squares  approach  to  minimize  a  functional,  S(x),  which  contains  two  terms  corre¬ 
sponding  to  the  two  set  of  observations  (N  receiver  function  sample  points,  y,  and 
M  dispersion  observations,  z)  adjusted  to  an  approximate  unit  variance.  A  control 
parameter,  p,  is  used  to  adjust  the  importance  of  each  set  of  observations. 

No  a  priori  approach  is  available  to  select  an  optimal  weighting  parameter,  and 
the  choice  is  made  best  after  a  systematic  study  of  the  inversions  produced  with 
a  range  of  weights  on  each  set  of  observations.  At  times  the  appropriate  value  of 
the  trade-off  parameter  can  be  made  by  studying  the  fits.  You  may  not  want  to  fit 
part  of  a  waveform  that  you  know  may  be  corrupted  by  off-azimuth  arrivals,  or  you 
may  not  want  to  fit  a  rough  portion  of  a  noisy  dispersion  cure.  These  seismological 
decisions  are  not  easily  coded,  but  the  decisions  are  easy  after  a  careful  study  of  the 
matches  to  the  signals  for  a  range  of  inversion  weights.  The  overriding  philosophy 
is  simply  to  construct  a  model  that  suitably  matches  the  observed  dispersion  and 
receiver  functions. 


Iterative  Jumping  Inversion  Scheme 


We  used  a  jumping  algorithm  to  jointly  invert  multiple  receiver  functions  (such  as 
those  with  different  bandwidths  or  corresponding  to  different  P-wave  incidence  an¬ 
gles)  and  surface-wave  observations  for  shear- wave  velocity  [Constable  et  ah,  1987]. 
The  jumping  scheme  allows  us  to  implement  a  smoothness  constraint  in  the  inver¬ 
sion  by  minimizing  a  model  roughness  norm  [Constable  et  al.,  1987;  Ammon  et  ah, 
1990]  that  can  trade-off  with  the  prediction  error.  Our  inversion  weighting  scheme 
accounts  for  the  different  number  of  data  points  and  different  physical  units  of  each 
data  set,  and  incorporates  an  a  priori  parameter  that  allows  to  investigate  the  rela¬ 
tive  influence  of  each  data  set  in  the  resulting  models  [Julia  et  al.  2000].  Our  model 
parameters  are  velocities  of  fixed-thickness  layers  overlying  a  half-space.  The  system 
of  equations  to  be  inverted  is 


(1.2) 
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\  a  A 
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where  Ds  and  Dr  are  the  partial  derivative  matrices  for  the  dispersion  measurements 
and  the  receiver  function  estimates,  respectively,  rs  and  rr  are  the  corresponding 
vectors  of  residuals,  m  is  the  vector  of  S-wave  velocities,  rn0  is  the  starting  model,  and 
A  is  a  matrix  that  constructs  a  depth-dependent  weighted  second  difference  of  the 
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model.  The  partial  derivative  matrices  and  the  vectors  of  residuals  are  normalized 
to  equalize  the  different  number  of  data  points  and  physical  units  in  the  data  sets. 
A  number  of  trade-off  parameters  must  be  specified  before  inversion.  The  influence 
factor,  p,  controls  the  trade-off  between  fitting  receiver  functions  and  dispersion 
curves,  and  the  smoothness  parameter,  a,  controls  the  trade-off  between  data  fitting 
and  model  smoothness.  The  parameter  q  =  1  -  p,  so  that  p  is  meaningless  outside 
the  range  0  <  p  <  1. 

We  begin  with  an  example  illustrating  the  inversion.  Figure  4  is  a  plot  of  the  ob¬ 
served  and  predicted  radial  receiver  functions  and  fundamental  mode  group  velocity 
values  for  station  SODA  located  near  the  western  edge  of  the  Arabian  Shield.  The 
receiver  function  sampled  the  structure  into  the  Shield  and  away  from  the  coast. 
The  resulting  shear-velocity  profile  is  shown  on  the  right.  The  shear-velocity  model 
layers  are  one-km  thick  for  depths  down  to  eight  km,  two-km  thick  from  8  to  50 
km  depth  and  five-km  thick  from  50  to  100  km  depth,  beneath  which  is  a  uniform 
half-space.  The  trade-off  parameters  for  the  inversion,  and  ,  were  chosen  empirically, 
after  a  systematic  study  of  the  inversion  results  produced  for  a  range  of  values,  as 
described  in  Julia  et  al.  [2000].  The  starting  model  was  a  uniform  8.0  km/s  P-wave 
velocity  half-space.  The  resulting  model  is  relatively  simple,  with  a  strong,  shallow 
velocity  gradient  with  a  two  layer  crust  underlain  by  a  transitional  lower  crust.  The 
upper  mantle  is  smooth,  but  contains  a  low  velocity  zone  below  50  km. 

Although  simple,  this  model  poses  some  interpretational  challenges.  First  the  shal¬ 
low  structure  appears  to  be  too  slow  when  compared  with  more  detailed  models  of 
shallow  shield  structure  which  suggest  an  average  velocity  of  3.2  km/s  for  the  top 
kilometer  [Mokhtar  et  al..  1988].  Our  surface  low-velocity  structure  is  3-to-4  km 
thick,  much  more  than  the  400  m  layer  in  higher-frequency  studies.  After  consider¬ 
able  numerical  experimentation  we  found  it  necessary  to  omit  the  shortest  periods 
for  Love  waves  (T  =  5  and  9  s)  to  remain  consistent  with  geologic  and  seismic 
constraints  on  shallow  structure  of  the  Shield  [e.g.  Mokhtar  et  ah,  1988]. 

Also  worth  close  inspection  is  the  need  for  the  decrease  in  velocity  throughout  the 
upper  mantle  in  the  model.  The  half-space  velocity  is  significantly  low  for  a  shield 
and  a  simple  check  of  the  long-period  dispersion  shows  values  that  are  much  lower 
than  those  commonly  observed.  The  models  deepest  velocities  are  controlled  by  the 
long-period  Rayleigh  waves,  and  not  surprisingly,  since  our  longest  periods  is  60  s, 
we  cannot  uniquely  resolve  much  of  the  deeper  structure  in  the  model.  However, 
although  our  observations  cannot  unambiguously  constrain  the  deep  structure,  they 
are  sensitive  to  deep  features  in  the  model,  particularly  an  infinitely  deep  half-space 
located  at  such  a  shallow  depth.  Our  concern  is  that  the  anomalous  structure  in 
the  mantle  may  still  trade  off  with  features  in  the  crust.  To  minimize  the  extent 


CHAPTER  1.  INTRODUCTION  &  METHOD  OVERVIEW 


15 


that  deep  velocity  features  may  influence  crustal  features,  vve  decided  to  constrain 
our  results  to  transition  smoothly  into  a  priori  mantle  structures.  Although  we  do 
not  know  a  prion  the  deep  structure,  we  are  confident  that  PREM  [Dziewonski 
and  Anderson,  1980]  or  slight  modifications  from  PRElVI  velocities  are  certainly  a 
better  choice  than  an  unconstrained  component  in  our  band-limited  inversion.  In 
practice  we  can  implement  constraints  to  be  consistent  with  global  aspherical  shear 
velocity  models  since  we  need  only  know  the  values  of  deeper  mantle  shear  velocity 
beneath  the  station  (which  we  can  get  for  the  shear-wave  models).  Of  course  we 
can  also  (and  plan  to)  include  longer-period  tomographic  dispersion  measurements 
from  other  studies  into  our  inversion. 


Including  a  priori  Velocity  Information 

The  upper  mantle.  The  upper  mantle  low  velocity  zone  in  Figure  5  may  be  a 
consequence  of  the  limited  resolving  power  at  depth  of  our  data  set.  To  test  the 
hypothesis  we  incorporated  a  priori  information  to  compel  the  deepest  layers  in  our 
model  to  resemble  values  from  global  mantle  models  like  S12WM13  model  [Su  et 
ah,  1994],  We  chose  to  constrain  our  solution  by  appending  the  following  set  of 
equations  to  the  original  system  (1)  [Jackson,  1972], 

(1.3)  W  771,  —  W  •  TTlapriorri 

where  W  is  a  diagonal  matrix  of  constraint  weights  and  m  contains  a  priori  prede¬ 
fined  velocity  values.  Equation  1.2  becomes 
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The  scalar  weight  e  is  redundant,  but  it  allows  us  to  adjust  the  importance  of  the  a 
pnori  velocity  weights  relative  to  the  model  fits  and  roughness  with  one  convenient 
parameter.  We  also  thicken  the  model  to  move  the  half-space  to  a  depth  of  500- 
km,  well  below  the  depth  of  influence  in  our  band  width.  The  resulting  models 
are  satisfactory  in  that  they  predict  dispersion  values  more  consistent  with  global 
observations  than  do  the  truncated-depth  models. 


Chapter  2 


Bulk  Properties  of  the  Continental 
Crust 


The  continental  crust  covers  approximately  41%  of  the  Earth’s  surface  and  accounts 
for  79%  of  the  volume  of  the  Earth’s  crust  [Condie,  1997].  Approximately  71%  of 
the  continental  crust  lies  above  the  mean  sea  level  [McLennan, 1992].  The  conti¬ 
nental  crust  is  a  composite  of  materials  fractionated  from  the  primordial  mantle 
“geochemical  reservoir”  by  geologic  processes  operating  during  the  4500  Ma  history 
of  the  Earth.  The  continental  crust  is  usually  divided  into  10  continents  and  four 
submerged  microcontinents.  The  mean  elevation  of  continental  crustal  material  is 
about  125  meters  above  sea-level  (including  the  29%  that  lies  beneath  the  ocean). 
The  continental  crust  ranges  in  thickness  from  about  10-80  km.  The  average  values 
for  the  crust  is  about  40  km  [Mooney,  1994],  although  prvious  estimates  of  35  or  38 
are  often  used  in  calculations.  The  crust  is  thickest  beneath  active  orogenie  regions, 
and  probably  the  thickest  crust  lies  beneath  Tibet  and  the  Bolivian  Altiplano.  In¬ 
cluding  the  ocean  sediment  of  continental  affinity,  the  volume  of  continental  crust 
is  estimated  to  be  approximately  7.35  x  109  km3  [McLennan,  1992].  The  lengthy 
history  of  the  crust  produces  a  complex  assembledge  of  materials  with  a  variety  of 
densities  and  elastic  parameters.  Assuming  a  mean  density  of  2800  kg/m3.  the  mass 
of  the  crust  is  approximately  2.06  x  1022  kg,  representing  about  0.35%  of  the  total 
mas  of  the  Earth  [McLennan,  1992]. 

The  single  outstanding  feature  of  continental  crustal  chemistry  is  the  enrichment 
of  the  incompatible  trace  elements.  From  the  concentrations  shown  in  Figure  2, 
we  can  also  see  that  the  amount  of  enrichment  correlates  well  with  the  affinity  of 
a  given  element  to  the  liquid  produced  by  a  mantle  melt.  The  bulk  chemistry  of 
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the  continental  is  approximately  that  of  tonalite,  with  approximately  65%  Si02 
[O’Nions,  1992].  In  general,  the  lower  crust  is  thought  to  be  more  mafic  than 
the  upper  crust,  indicating  an  increase  importance  of  basalt  in  the  lower  crust. 
Crustal  xenoliths  strongly  indicate  that  the  lower  crust  is  more  mafic  than  the 
upper  crust.  Exposed  sections  of  the  lower  crust  agree  with  this  general  trend, 
but  are  not  as  basaltic  in  character  as  the  samples  found  in  xenoliths  [Rudnick, 
1992].  The  continental  crust  is  the  only  repository  of  information  for  a  major  part 
of  the  earth  history  [McLennan,  1992]  since  the  relatively  young  age  of  the  oceans 
(200  Ma)  supplies  information  on  earth  evolution  for  only  a  small  fraction  of  the 
age  of  the  Earth.  Models  of  crustal  evolution  are  primarily  based  on  geochemical 
(isotopic)  evidence  although  recent  efforts  at  global  surveys  using  seismic  data  are 
adding  new  information  to  the  problem  [Nelson,  1991;  Durrheim  and  Mooney,  1991 
and  1994].  Early  in  Earth  history,  the  planet  had  very  little,  if  any,  crust.  Continued 
bombardment  of  by  planetismals  and  debris  remaining  from  solar  system  formation 
continually  broke  up  the  slag  forming  above  the  mantle.  Eventually,  the  rate  of 
bombardment  slowed  enough  to  allow  the  formation  of  crust.  We  have  samples  of 
continental  crustal  material  back  to  about  3800  Ma.,  so  only  the  first  700  Ma  are 
without  record  [O’Nions,  1992]  (the  oldest  terrestial  material,  zircons  from  western 
Australia  have  been  dated  to  approximately  4200  Ma.).  The  estimated  mean  age  of 
the  crust  is  about  2200  300  Ma. 


Poisson’s  Ratio 

Poisson’s  Ratio,  a,  is  an  elastic  parameter  that  is  often  used  to  characterize  different 
rock  types.  Physically,  Poisson’s  Ratio  can  be  defined  if  you  consider  an  experiment 
stretching  an  elastic  rod  of  diameter  d  and  length  L.  As  you  stretch  the  rod,  the 
length  and  the  diameter  will  change.  Poisson’s  ratio  is  then  defined  as 

A  d/d 

°  =  a  ~l]l' 

For  most  materials  as  you  stretch  the  rod  the  diameter  decreases,  and  so  o  is  positive. 
For  most  rocks,  typical  values  of  Poisson’s  Ratio  is  on  the  order  of  0.25,  so  the 
fractional  change  in  diameter  is  about  1/4  that  of  the  fractional  change  in  length 
when  the  material  is  stretched. 

For  seismic  applications,  Poisson’s  Ratio  can  be  related  to  the  seismic  velocities  of 
P  and  S  waves,  which  can  be  measured.  Let  Vp  represent  the  P-velocity  and  Vs 
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represent  the  S  velocity.  Then 


Vp 

Vs 


2(1  -a) 

1  -  2a 
1 

1  —  2a  " 


1. 


The  inverse  relationship  is 


a 


The  trick  then  is  to  relate  a  to  rock  type.  The  relationship  between  a  and  rock 
composition  is  not  unique  by  any  metric,  but  the  information  can  be  used  in  some 
instances  to  distinguish  between  different  rock  types.  Even  better  is  knowledge  of 
both  Vp,  and  Vs.  The  relationship  between  Vp/Vs  and  Poisson’s  ratio  is  one-to-one, 
so  the  two  can  be  used  interchangeably. 


Receiver  Functions  for  Simple  Earth  Models 


In  relatively  simple  structures,  the  waveform  is  a  composite  of  P-to-S  converted 
waves  that  reverberate  in  the  structure  beneath  the  seismometer.  Modeling  the  am¬ 
plitude  and  timing  of  those  reverberating  waves  can  supply  valuable  constraints  on 
the  underlying  geology.  If  the  main  structural  boundaries  are  approximately  sub¬ 
horizontal  in  the  region  of  the  recording  station,  the  main  features  of  the  structure 
can  be  approximated  using  a  sequence  of  nearly-horizontal  interfaces  bounding  re¬ 
gions  with  strong  contrasts  in  elastic  properties.  In  that  case,  the  arrivals  generated 
by  each  sharp  (that  is,  sharp  relative  to  the  shortest  wavelength  in  the  observations) 
looks  something  like  the  signal  shown  in  Figure  2.1. 

In  more  complex  structures,  the  arrivals  may  also  include  scattered  surface  waves, 
and  the  response  may  be  complicated. 

The  travel  times  of  these  waves  can  be  simply  expressed  in  terms  of  the  horizontal 
and  vertical  slownesses,  p  and  r)a  and  The  horizontal  slowness,  p  is  also  called 
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Direct  P  wave 


Figure  2.1:  Canonical  receiver  function  for  layer  over  a  half-space  model. 


the  ray  parameter,  and  can  be  estimated  from  the  slope  of  travel-time  curves.  The 
vertical  slownesses  are  related  to  the  seismic  velocity  and  the  ray-parameter  by 


V0 


1 


Vs2 


and 


Pa  = 


Vp2 


P 


In  a  receiver  function,  as  a  consequence  of  the  deconvolution,  all  times  are  specified 
relative  to  the  initial  P-wave  arrival.  Then 


Tps 

TppPmS 


Ti 


PsPmS 


=  T 


PpSPmS 


{Pd  ~  Pa)  ■  h 

(f]0  +  Va)-  h 

2  •  P0  ■  h 


The  amplitudes  of  the  arrivals  in  a  receiver  function  depend  on  the  incidence  an¬ 
gle  of  the  impinging  P-wave  and  the  size  of  the  velocity  contrasts  generating  the 
conversions  (Ps)  and  multiples  (PsPhs,  PpShs,  PpPhs)  .  The  arrival  times  of  the 
converted  phase  and  multiples  depend  on  the  depth  of  the  velocity  contrast,  the 
P  and  S  velocity  between  the  contrast  and  the  surface,  and  the  P-wave  incidence 
angle,  or  ray  parameter.  The  relative  amplitudes  of  the  later  arrivals  and  their  fre¬ 
quency  content  depend  on  the  nature  of  the  velocity  transition  -  e.g.  is  it  abroad  or 
relatively  sharp  velocity  change? 
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PpSms 


Station 
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Figure  2.2:  Ray  paths  for  the  main  arrivals  in  a  receiver  function  generated  by  a 
simplified  crustal  model.  Each  cross  section  is  cut  along  the  plane  containing  the 
receiver  and  the  source.  The  left  panel  shows  all  the  rays  for  the  largest  arrivals, 
the  panels  on  the  right  separate  the  converted  phases  from  the  multiples  for  clarity, 
for  clarity. 
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Poisson-Ratio  Estimation  Using  Amplitude  Stacking 

This  Poisson’s  Ratio  panel  in  Lithospheric  Structure  in  Africa  and  Western  Eura¬ 
sians  a  simple  implementation  of  the  h  —  k  stacking  algorithm  of  Zhu  and  Kanamori 
[2000].  h  represents  thickness,  and  k  represents  the  Vp/Vs  ratio.  A  grid-search 
is  performed  for  a  range  of  thickness  and  Poisson  Ratio’s  assuming  a  constant  P- 
velocity.  For  each  thickness  and  velocity,  the  arrival  times  of  the  Ps,  PpPrnP , 
PsPmS ,  and  PpSrnS  phases  (Figure  2.1)  are  computed  and  the  receiver  function 
amplitudes  corresponding  to  these  times  are  summed.  The  optimal  estimates  for  the 
values  of  these  parameters  is  that  which  produces  the  maximum  amplitude.  Since 
the  second  set  of  multiples  PsPmS ,  and  PpSmS  have  the  opposite  polarity  of  the 
Ps  and  PpPmS  multiple,  the  former  amplitude  is  subtracted  to  the  others.  Also, 
since  the  quality  of  the  multiples  varies  and  is  generally  not  as  good  as  the  Ps  phase 
quality,  a  weighted  sum  of  these  amplitudes  is  employed  If  we  let  P(h,a)  represent 
the  sum  for  a  particular  thickness,  h>  and  Poisson’s  Ratio,  a ,  then  we  have 

(2.1)  P(h,(j)  =  WpsR(Tps)  +  W pppmS H(TppPms)  ~  W PpSrnS H{TppSms) 

where  the  Wfs  represent  weights  and  R(t )  represents  the  receiver  function.  A  car- 
toonish  vision  of  the  ray  paths  are  shown  in  Figure  2.2.  Note  that  implicit  in  the 
assumptions  of  this  method  is  a  simple,  uniform  crust  underlain  by  a  flat  crust- 
mantle  transition. 

The  choice  of  weights  in  Equation  (2.1)  is  subjective.  Most  values  are  adjusted 
by  trial  and  error  with  an  eye  on  sensitivity.  We  use  a  slightly  different  approach 
(Figure  2.3)  in  that  we  examine  the  stacking  images  for  each  of  the  three  phases 
used  to  produce  a  common  receiver  function  h  —  k.  stack.  With  this  approach,  we 
can  determine  the  consistency  of  the  phases  used  in  the  stacks.  We  also  examine  the 
moveout  of  the  different  phases  as  a  function  of  ray  parameter.  If  one  of  the  phases 
corresponds  to  a  poorly  defined  feature,  we  can  decrease  it’s  weight.  For  example, 
in  Figure  2.3,  the  Ps  and  PpPmS  phases  are  consistent,  but  the  later  multiple 
arrivals  are  inconsistent  with  the  other  two.  This  suggests  not  using  the  second  set 
of  multiples  ( PsPmS  and  PpSmS)  in  the  stacking.  We  have  run  bootstrap  error 
estimates  on  many  of  the  data,  and  together  with  the  consistency  of  the  results  and 
a  comparison  based  with  all  of  the  stations  that  we  have  analyzed,  we  grade  each 
estimate  across  the  range  of  A,  B,  C,  D,  and  F.  While  inherently  fuzzy,  these  grades 
provide  experienced  quality  assessments  for  those  wishing  to  use  the  results. 

To  perform  the  stacking  we  must  also  estimate  the  vertical  slownesses,  which  depend 
on  the  ray  parameter  (computed  from  travel-time  curves),  and  the  assumed  average 
P-wave  speed  for  the  crust.  The  results  are  sensitive  to  the  assumed  velocity  so  the 
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Figure  2.3:  Sample  receiver-function  stacking  example. 
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speed  that  you  use  is  important.  We  used  the  average  speed  obtained  at  each  station 
from  Crust  2.0  (actually  a  slowness  average),  which  allows  a  direct  comparison 
between  our  results  and  the  reference  model,  and  keeps  our  results  as  consistent 
with  the  global  reference  as  possible. 


Poisson’s  Ratio  &;  Crustal  Thickness  Estimates 


In  this  section  we  present  a  summary  view  of  the  results  of  receiver-function  stacking 
to  estimate  the  crustal  thickness  and  Poisson’s  ratio  across  the  study  region.  All  of 
our  results  are  assigned  a  quality  rank,  which  is  listed  along  with  all  estimates  in  ap¬ 
pendix.  Here,  we  summarize  the  results  graphically  and  compare  our  measurements 
with  those  from  previous  studies. 


Comparison  With  Zandt  and  Ammon  [1995].  Values  for  crustal  thickness 
and  Poisson’s  Ratio  estimated  by  Zandt  and  Ammon  [1995]  are  summarized  in 
Figure  2.4.  The  study  includes  76  measurements,  and  was  performed  by  picking 
the  arrival  times  of  the  converted  and  multiple  phases,  not  using  the  less  subjec¬ 
tive  stacking  approach  developed  later.  Experience  suggests  that  careful  picking 
produces  results  compatible  with  the  stacking  results.  Zandt  and  Ammon  sepa¬ 
rated  their  results  into  crustal  provinces  guided  by  the  map  of  basement  age  found 
in  Condie  [1993].  The  results  for  most  provinces  included  much  scatter  with  the 
exception  of  shields,  which  were  uniformly  high  in  Poisson’s  ratio.  The  average 
thicknesses  for  each  crustal  age  province  produced  a  systematic  variation  with  age 
as  shown  in  Figure  2.4.  The  results  of  the  study  were  based  predominately  on  two 
large  earthquakes  from  1994  (the  large,  deep  Bolivia  and  Fiji  events)  and  were  heav¬ 
ily  biased  to  North  America  -  a  result  of  the  global  seismic  network  configurations 
at  the  time.  Only  a  handful  of  stations  are  common  to  the  present  analysis  and  the 
1995  study,  so  direct  quantitative  comparisons  are  unavailable. 


Comparisons  With  Crust  2.0.  One  of  the  more  interesting  comparisons  that 
we  can  make  is  with  the  predominantly  refraction  (P-wave  speed)  constructed  model 
Crust  2.0  [Bassin  et  aL,  2000].  Crust  2.0  is  a  model  widely  used  in  global  and  regional 
seismic  waveform  modeling,  and  forms  the  basis  for  crustal  corrections  in  long-period 
shear-velocity  imaging  efforts.  Crust  2.0  contains  a  detailed  map  of  crustal  thickness, 
but  has  a  low  spatial  resolution  of  variations  in  Poisson’s  Ratio  or  equivalently,  the 
Vp/Vs  ratio.  Thicknesses  are  compared  in  Figure  2.6,  which  is  useful  only  for  first- 
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Figure  2.4:  Global  variations  in  Poisson’s  ratio  observed  b}^  Zandt  and  Ammon 
(1995).  The  simplified  geologic  map  is  based  on  REFERENCE. 


order  comparison  of  the  two  sets  of  results.  Many  observations  agree  to  within  2.5 
km,  probably  the  resolution  of  either  method,  and  most  of  the  measurements  are 
within  about  3.5  kilometers  of  the  values  listed  in  Crust  2.0.  About  90%  of  the 
observations  are  within  7  kilometers.  As  in  Zandt  and  Ammon  [1995],  observations 
are  separated  into  crustal  province  to  illuminate  age-dependent  differences  (which 
are  not  readily  clear).  The  thickness  values  correlate  reasonably  well  (Figures  2.5 
and  Figures  2.6  ).  Our  results  for  the  Tanzania  Craton  suggest  a  thinner  crust  than 
that  contained  in  Crust  2.0,  but  the  other  values  are  more  consistent  such  as  those 
in  the  Arabian  Shield.  Thicker  crusts  are  observed  in  the  shields  in  west-central 
Africa  and  on  the  western  edge  of  the  Tanzania  Craton.  The  thinnest  crusts  are 
observed  across  north-central  Europe  (less  than  30  km),  the  thickest  in  mountain 
belts,  including  the  Urals,  Tian  Shan,  and  Pamir  region. 

In  Fig  ure  2.7  we  compare  the  our  Vp/Vs  ratio  estimates  with  those  included  in  Crust 
2.0.  Crust  2.0  has  an  extremely  simple  variation  in  Poisson’s  ratio  (there  are  few 
non  receiver-function  based  estimates  of  this  quantity  [Zandt  and  Ammon,  1995]. 
Individual  results  are  uncorrelated,  but  are  bounded  by  the  range  of  values  assumed 
in  Crust  2.0.  Table  2.1  summarizes  our  results  by  crustal  age  province,  and  we 
include  the  Crust  2.0  values  sampled  at  the  same  locals  for  reference.  Quantitative 
comparison  between  the  two  populations  is  difficult,  the  information  is  simply  not 
contained  in  Crust  2.0. 

For  completeness,  we  also  summarize  the  results  geographically,  using  the  maps 
shown  in  Figure  (2.8).  The  sparsity  of  stations  makes  drawing  conclusions  difficult. 
But,  as  mentioned  earlier,  the  thinnest  crusts  are  located  in  north-central  Europe, 
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Table  2. 


Crustal  Thickness  &  Poisson’s  Ratio  By  Age  Province. 
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Figure  2.5:  Observed  differences  between  the  crustal  thickness  estimated  using  re¬ 
ceiver  function  stacking  and  those  estimated  in  the  construction  of  model  Crust  2.0. 
The  diagonal  lines  indicate  on  standard  deviation  in  the  difference  between  the  two 
crustal  thickness  estimates.  Horizontal  trends  usually  indicate  stations  within  the 
same  cells  in  model  Crust  2.0. 
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Crustal  thickness  in  Crust2.0  -  Crustal  Thickness  in  This  Study 

Figure  2.6:  Distribution  of  differences  between  the  crustal  thickness  estimated  using 
receiver  function  stacking  and  those  estimated  in  the  construction  of  model  Crust 
2.0.  The  average  shown  in  the  chart  is  not  adjusted  for  crustal  province,  so  it  only 
useful  as  a  comparison  of  the  results  samples  as  the  crust  is  in  this  study. 
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Figure  2.7:  Observed  of  differences  between  the  the  Vp/Vs  ratio  estimated  using 
receiver  function  stacking  and  those  contained  in  model  Crust  2.0. 


the  thickest  beneath  active  orogens.  Shield  values  are  intermediate,  but  slightly 
thicker  in  western  Africa.  A  similar  summary  for  Vp/Vs  ratios  is  presented  in 
Figure  2.9. 

We  also  performed  a  literature  search  of  results  from  across  the  globe,  including 
receiver-function  estimates  of  Vp/Vs  ratios  from  North  America,  Japan,  Australia. 
The  results  are  shown  in  Figure  2.10,  which  also  contains  a  summary  of  studies 
already  discused.  The  histogram  labeled  Holbrook  et  al.  (1992)  include  thickness 
values  from  about  90  active-source  seismic  experiment,  but  included  only  11  Vp/Vs 
ratio  estimates.  One  remarkable  result  is  the  systematic  crustal  thickness  pattern  ob¬ 
served  from  the  different  receiver  function  observations  -  crustal  thickness  increases 
by  2-5  km  with  the  addition  of  sediments  on  top  of  shield  to  form  Platforms,  then 
decreases  beneath  the  Paleozoic  orogens.  The  patterns  beneath  the  active  Mesozoic- 
Cenozoic  regions  are  much  more  difficult  to  simplify,  as  the  scatter  of  values  within 
that  provinces  is  large.  The  unusual  nature  of  the  Tibetan  crust  is  highlighted  in 
the  global  summary. 
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Figure  2.8:  Observed  crustal  thicknesses  from  receiver  function  stacking. 
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Figure  2.9:  Observed  crustal  Vp/Vs  ratios  from  receiver  function  stacking. 
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Figure  2.10:  Graphical  comparison  of  results  from  earlier  studies,  this  work,  and  a 
search  of  the  scientific  literature  for  receiver  function  stacking  estimates  across  the 
globe. 


Chapter  3 


Shear- Velocity  Estimation 


As  described  in  the  introduction,  subsurface  geology  generally  has  a.  broad  wavenum¬ 
ber  spectrum  containing  spatially  localized  broadband-wavenumber  changes  in  ve¬ 
locity  near  Earths  major  geologic  boundaries  and  smooth  low-wavenumber  variations 
in  regions  of  relatively  uniform  geologic  structure.  Access  to  the  full  spectrum  of 
earth  structure  requires  that  we  exploit  signals  that  span  a  wide  frequency  range 
and  that  are  sensitive  to  the  entire  spectrum  of  heterogeneity.  Surface- waves,  travel 
times,  and  direct-wave  amplitudes,  for  example,  are  sensitive  to  smooth  variations 
in  earth  structure;  reflected  and  converted  waves  are  sensitive  to  velocity  contrasts. 
Combining  seismic  data  in  joint  inversions  is  an  obvious  approach  to  improve  es¬ 
timates  of  earth  structure.  To  successfully  combine  data  in  an  inversion,  we  must 
insure  that  all  the  data  are  sensitive  to  the  same  (or  related)  physical  quantities  and 
that  they  sample  or  average  structure  over  comparable  length  scales.  Advances  in 
surface-wave  tomography  have  provided  an  opportunity  to  combine  localized  surface- 
wave  dispersion  estimates  with  other  data  such  as  P-  and  S-wave  receiver  functions. 

Surface- wave  dispersion  measurements  are  sensitive  to  broad  averages,  or  low  wave- 
number  components  of  earth  structure.  They  provide  valuable  information  on 
the  absolute  seismic  shear  velocity  but  are  relatively  insensitive  to  sharp,  high- 
wavenumber  velocity  changes.  Generally  surface- wave  inversions  must  be  constrained 
using  a  particular  layer  parameterization  (e.g.  near-surface,  upper-crust,  lower  crust, 
mantle  lid,  deep  mantle),  resemble  an  a  priori  model,  or  substantially  smoothed 
to  stabilize  earth-structure  estimation.  Despite  these  drawbacks,  surface-wave  dis¬ 
persion  values  contain  important  constraints  on  the  subsurface  structure,  and  the 
general  increase  in  depth  sensitivity  with  depth  allows  an  intuitive  understanding  of 
their  constraints  on  structure.  Additionally,  modeling  dispersion  values  facilitates  a 
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broadband  inversion  by  reducing  the  dominance  of  Airy  phases,  which  pose  prob¬ 
lems  when  constructing  broad-band  misfit  norms  to  model  seismograms  directly. 
Perhaps  most  important  for  our  application  is  the  ability  to  localize  Earths  dis¬ 
persion  properties  using  seismic  tomography.  The  idea  is  now  well  established  and 
global  dispersion  models  exist  for  a  broad  range  of  frequencies  (e.g.  Larson  and 
Ekstrrn,  2001;  Stevens  et  ah,  2001).  The  localization  of  dispersion  allows  us  to 
isolate  the  variations  in  properties  spatially  and  global  models  of  surface-wave  dis¬ 
persion  exist  and  are  readily  available  for  application  to  other  studies  such  as  the 
proposed  work.  We  employ  two  dispersion  models  in  our  inversions  -  a  low-frequency 
model  of  Harvard  (Larson  and  Ekstrom,  2001)  and  the  Lawrence  Livermore  Model 
of  Michael  Pasyanos  (Pasyanos  and  Walter,  2001.  personal  communications  -  for 
updates).  The  Harvard  models  contain  smooth  variations  as  a  function  of  period, 
but  have  no  information  below  periods  of  30  seconds  (important  for  upper  crustal 
constraints).  The  LLNL  model  has  shorter  period  information,  but  can  be  rough 
from  period-to-period.  We  generally  apply  a  simple  smoothing  operator  to  the  raw 
curves  before  using  the  LLNL  results  in  an  inversion. 

Receiver  function  waveforms  are  a  composite  of  P-to-S  (or  S-to-P)  converted  waves 
that  reverberate  within  the  structure  near  the  seismometer.  Modeling  the  ampli¬ 
tude  and  timing  of  those  reverberating  waves  can  supply  valuable  constraints  on  the 
underlying  geology.  In  general,  the  receiver  functions  sample  the  structure  over  a 
range  of  10s  of  kilometers  from  the  station  in  the  direction  of  wave  approach  (the 
specific  sample  width  depends  on  the  depth  of  the  deepest  contrast).  Stations  sited 
near  geologic  boundaries  can  produce  different  responses  for  different  directions. 
Recent  innovations  in  receiver  function  analysis  include  more  detailed  modeling  of 
receiver  function  arrivals  from  sedimentary  basin  structures  (e.g.  Clitheroe  et  al., 
2000),  anisotropic  structures  (e.g.  Levin  and  Park,  1997;  Savage  1998),  estima¬ 
tion  of  Poissons  ratio  (e.g.  Zandt  et  al.,  1995;  Zandt  and  Ammon,  1995;  Zhu  and 
Kanamori,  2000,  Ligorra,  2000),  reflection-like  processing  of  array  receiver  functions 
(e.g.  Chevrot  and  Girardin,  2000;  Ryberg  and  Weber,  2000)  and  joint  inversions 
(e.g.  Ozalaybey  et  al.,  1997;  Du  and  Foulger,  1999;  Julia  et  al.,  2000). 

Ammon  and  Zandt  (1993)  used  surface- wave  dispersion  observations  to  try  and 
distinguish  between  competing  models  of  the  Mojave  desert,  but  Ozalaybey  et  al., 
(1997)  pioneered  a  formal,  joint  inversion  of  these  data.  They  nicely  illustrated  the 
value  of  even  a  limited  band  of  dispersion  values  to  help  reduce  the  trade-off  between 
crustal  thickness  and  velocity  inherent  in  receiver  function  analyses.  Specifically, 
they  used  Rayleigh-wave  phase  velocities  in  the  2025  second  period  range  to  help 
produce  stable  estimates  of  crustal  thickness  in  the  northern  and  central  Basin  & 
Range.  The  limited  bandwidth  did  not  permit  resolution  of  details  in  the  crust 
and  they  limited  their  inversion  (or  at  least  their  interpretation)  to  depths  above  40 
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km.  More  recent  authors  have  exercised  the  approach  and  combined  the  data  with 
additional  a  priori  model  constraints  (Du  and  Foulger,  1999;  Julia  et  ah,  2000). 
Recent  accomplishments  in  global  and  regional  tomography  now  provide  a  more 
complete  band  of  dispersion  measurements  to  combine  with  receiver  functions  that 
allow  us  to  improve  the  resolution  of  earlier  works. 

Seismic  inversion  is  not  a  simple  exercise.  Each  result  requires  performing  dozens 
of  inversions  to  assess  the  stability  of  the  result  and/or  identify  the  robust  features 
in  the  solution.  For  the  inversion  of  dispersion  curves,  we  often  experimented  with 
adjusting  the  bandwidth  of  the  observations  or  excluding  Love  or  Rayleigh  waves  to 
explore  the  consistency  and  sensitivity  of  those  data.  After  months  of  experiments 
and  exploration  of  these  issues  we  found  that  most  of  the  time  it  was  prudent  to 
often  ignore  the  Love  wave  measurements,  which  are  less  well  constrained  in  the 
tomographic  models  than  the  corresponding  Rayleigh  wave  dispersion  values.  The 
extent  of  the  differences  between  these  values  is  an  interesting  scientific  target,  but 
we  do  not  believe  that  simply  extracting  curves  from  the  models  and  showing  that 
it  is  difficult  to  fit  them  simultaneously  is  sufficient  documentation  of  the  prevalence 
of  anisotropy.  The  tomography  for  Love  and  Rayleigh  are  decoupled  and  the  cov¬ 
erage  and  resolution  of  each  differs.  The  result  may  be  artificial  differences  in  the 
dispersion  values.  What  we  can  state  is  that  we  had  a  difficult  time  fitting  Love  and 
Rayleigh  dispersion  curves  simultaneously.  But  we  could  fit  the  receiver  function 
with  either  curve. 


An  Example  Inversion  In  the  joint  inversion  of  a  waveform  with  a  dispersion 
curve,  we  also  have  to  deal  with  the  weighting  of  the  two  different  data  sets.  If  the 
data  constrain  orthogonal  components  of  the  model,  the  results  are  not  too  sensitive 
to  the  relative  weights  applied  to  the  data  sets.  We  have  taken  two  approaches  to 
solving  this  problem,  first  we  have  run  a  large  number  of  inversions  with  different 
weighting  factors  to  identify  an  influence  parameter  (relative  weight)  for  the  two  data 
sets.  Second,  we  have  used  a  multiple  objective  optimization  approach  [Ammon, 
2003],  to  map  out  a  more  objective  weight  selection  approach.  Both  analyses  suggest 
roughly  the  same  value.  For  ease  of  reference,  we  repeat  equation  3  here 
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where  the  subscript  r  identifies  receiver  function  quantities  and  the  subscript  s 
identifies  surface-wave  quantities,  a  and  e  control  the  importance  of  smoothing  and 
a  prion  model  constraints. 
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Figure  3.1:  Location  map  for  station  TAM,  used  in  the  example  inversions. 


As  an  example,  we  present  the  several  inversions  of  the  data  recorded  at  station 
TAM,  Tamanrassett,  Algeria  (Figure  3.1).  TAM  is  located  within  the  Hoggar  swell, 
an  uplifted  region  in  north-central  Africa.  Observations  from  the  southwest  are 
particularly  clear  at  TAM  and  we  use  those  for  our  illustration.  We  used  a  simple 
crustal  model  resting  on  a  PREM  mantle  as  the  initial  model,  and  ran  each  inver¬ 
sion  for  5  iterations.  Although  we  applied  smoothing  in  the  inversion,  we  did  not 
incorporate  any  a  prion  information  other  than  requiring  the  model  to  blend  in  to 
PREM  below  a  depth  of  200  km.  Figure  3.2  shows  the  results  for  a  receiver-function 
dominated  inversion  (influence  parameter  p  =  0.05,  q  =  1  —  p  =  0.95).  Love  wave 
dispersion  values  are  shown  for  reference  but  were  not  used  in  the  inversion.  Also 
shown  is  a  low-frequency  receiver  function,  which  also  was  left  out  to  simplify  the 
inversion  (but  we  have  the  capability  and  usually  include  different  bandwidth  re¬ 
ceiver  functions  in  the  inversions).  Even  with  most  of  the  weight  on  the  receiver 
function,  the  fit  to  the  Rayleigh  waves  is  reasonable.  The  solution  obtained  after 
reversing  the  weights  is  shown  in  Figure  (3.3).  The  result  here  is  smoother  and  the 
fit  to  the  Rayleigh-wave  dispersion  values  is  excellent,  but  the  fit  to  the  receiver 
function,  which  requires  some  sharp  boundaries,  has  degraded. 

We  can  perform  a  systematic  variation  in  the  influence  parameter  to  compute  an 
optimal  weight.  The  results  of  the  systematic  search  are  shown  in  Figure  3.4.  The 
chart  is  a  standard  multiple  objective  optimization  misfit  diagram,  that  shows  the 
ideal  misfits  and  the  misfits  found  for  a  range  of  influence  parameter  (labeled  points). 
The  point  closest  to  the  ideal  fit  is  roughly  0.40.  Inversion  results  for  this  balanced 
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Influence  Parameter  =  0.05 


Period  (s) 

Figure  3.2:  Sample  inversion  showing  the  fits  to  a  broadband  receiver  function  and 
Rayleigh  dispersion.  In  this  example  the  influence  parameter  is  0.05,  which  produces 
an  inversion  dominated  by  the  receiver  function.  The  model  outlined  in  green  is  that 
obtained  from  a  more  balanced  weight. 
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Influence  Parameter  =  0.95 


Period  (s) 

Figure  3.3:  Sample  inversion  showing  the  fits  to  a  broadband  receiver  function  and 
Rayleigh  dispersion.  In  this  example  the  influence  parameter  is  0.95,  which  produces 
an  inversion  dominated  by  the  dispersion  (in  this  case,  the  Rayleigh  waves).  The 
model  outlined  in  green  is  that  obtained  from  a  more  balanced  weight. 
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Receiver  Function  Misfit 


Figure  3.4:  Mulitple  objective  optimization  misfit  diagram  showing  the  variation 
in  misfit  to  the  two  data  sets  as  for  a  range  of  influence  parameters.  The  optimal 
influence  parameter  produces  misfits  closest  to  an  ideal  solution  fit,  that  is  identified 
by  the  the  point  with  coordinates  given  by  the  lowest  misfit  for  each  data  set. 
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Influence  Parameter  =  0.4 


Period  (s) 

Figure  3.5:  Sample  inversion  showing  the  fits  to  a  broadband  receiver  function  and 
Rayleigh  dispersion.  In  this  example  the  influence  parameter  is  0.40,  which  produces 
a  balanced  inversion. 
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weight  are  shown  in  Figure  3.5.  The  fit  to  both  data  sets  in  this  ease  is  superb 
(probably  better  than  warranted)  and  the  model  is  relatively  smooth,  but  contains 
sharp  features  where  needed,  such  as  near  the  surface  and  near  the  crust-mantle 
boundary. 


Shear- Velocity  Models 


Earth  is  not  always  amendable  to  simplified  models  of  its  structure.  We  have  pro¬ 
duced  a  library  of  one-dimensional  shear  velocity  inversions  to  the  lithosphere.  These 
models  are  at  best  approximations  to  the  true  structure  and  should  be  thought  of 
as  only  a  first  approximation  to  what’s  really  at  depth.  But  computer  programs 
always  provide  and  answer,  but  the  answers  vary  in  reliability.  In  an  appendix  we 
provide  the  results  of  all  171  shear-wave  speed  inversions.  Not  all  of  the  results 
are  reliable.  We  have  ranked  the  receiver  function  quality  in  the  tables  contained 
in  the  appendices.  Before  using  any  results  you  should  check  the  quality  of  the 
data.  Checking  the  quality  of  the  dispersion  values  is  more  difficult,  since  they  are 
provided  by  another  group  and  the  inherent  difficulties  associated  with  estimating 
reliable  uncertainties  in  tomography.  These  issue  complicate  the  use  of  our  shear- 
velocity  models,  and  necessitate  that  people  using  the  information  must  have  a  basic 
understanding  of  the  nature  of  geophysical  inverse  problems  and  some  understanding 
of  seismic  wave  propagation. 

The  shear-velocity  models  can  depend  strongly  on  the  group  velocities  used  in  the 
inversion.  We  have  performed  inversions  using  both  the  LLNL  [Pasyanos  and  Walter, 
2002]  and  the  Harvard  models  Larson  and  Ekstrom,  2001].  When  the  dispersion 
models  differ  substantially,  our  estimates  of  the  absolute  velocities  will  diverge; 
velocity  contrasts,  controlled  by  the  receiver  functions,  will  be  similar,  but  their 
depth  may  change  to  maintain  the  time  lags  of  the  corresponding  arrivals  observed 
on  the  receiver  functions.  Convergence  in  tomographic  models  to  uniform  dispersion 
models  is  the  only  way  to  remove  such  ambiguity. 


Relatively  Simple  Radial  Receiver  Function  Observations.  A  number  of 
stations,  predominantly  on  shields,  exhibit  relatively  simple  receiver  function  re¬ 
sponses  that  are  dominated  by  relatively  sharp  crust-mantle  boundary  arrivals.  We 
list  the  simplest  stations  in  Table  3.1  since  these  may  provide  ideal  stations  for  stud¬ 
ies  needing  to  observe  clear  teleseismic  observations  for  magnitude  or  for  teleseismic 
depth  estimation.  One  must  remember  that  almost  all  receiver  responses  depend  to 
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some  extent  on  the  incidence  and  angle  and  back  azimuth,  so  some  care  will  always 
be  needed  using  these  observations.  Note  that  networks  beginning  with  the  letter 
X  were  temporary  -  the  IRIS  PASSCAL  deployments. 


Table  3.1:  Stations  With  Simple  Radial  Receiver  Function  Observations 


Network  Code 

Station  Code(s) 

AS 

KAAO 

cz 

DPC 

G 

ECH,  SSB,  TAM 

GE 

MALT,  MTE,  STU, 

GT 

BGCA 

II 

ESK.  RAYN, 

IU 

KMBO,  KOWA,  MSKU 

KZ 

CHK,  ZRN 

MN 

VSL,  YVDD  W 

XD  (Tanzania) 

AMBA,  KOND,  LONG,  MB  WE,  MITU, 
MTOR,  PUGE,  RUNG 

XI  (Arabia  Shield) 

AFIF,  HALM,  RANI,  SODA, 

Exceedingly  Complex  Receiver  Function  Observations.  We  list  stations 
with  exceedingly  complex  receiver  functions  in  Table  3.2.  Models  derived  from 
these  data  are  unreliable  and  should  not  be  used.  In  a  few  instances  the  complexity 
may  be  the  result  of  sparse  data,  but  some  are  simply  indicative  of  seismic  stations 
located  in  complex  environments,  and  several  are  island  stations,  where  the  structure 
often  has  strong  lateral  variations.  Other  reasons  for  complexity  may  be  proximity 
to  large  basins  such  as  that  found  near  station  SFUC  in  southwestern  Spain. 


Table  3.2:  Stations  With  Exceedingly  Complex  Structure 


Network  Code 

Station  Code(s) 

GE 

APEZ,  MELI.  SFS,  SFUC,  ZKR 

MN 

AIO,  AQU,  CEL,  CUC,  GFA,  IDI,  TUE,  VAE 

XG 

HAMR,  KARS,  TTB01 

XI 

RIYD 
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Near-Surface  Complexity  Observations.  Receiver  functions  are  also  quite  use¬ 
ful  for  identifying  regions  with  near-surface  structural  complexity.  Most  of  the  com¬ 
plex  stations  fall  into  this  category,  since  they  invariably  have  long-duration  radial 
(and  tangential)  receiver  functions  that  ring  on  for  tens  of  seconds.  Some  observa¬ 
tions  also  show  substantial  near  surface  complexity  that  lasts  less  than  10  seconds. 
This  is  generally  an  indication  of  near-receiver  reverberations  in  a  shallow  structure 
underlain  (or  containing)  strong  elastic  contrasts.  At  times  the  indication  is  simply 
a  shift  in  the  P-wave  arrival  in  the  receiver  function,  other  times  the  P-wave  is  over- 
whelmened  by  a  very  strong  Ps  conversion  from  a  less  than  a  few  kilometers  depth. 
Stations  whose  receiver  functions  exhibit  these  characteristics  are  list  in  Table  3.3. 


Table  3.3:  Stations  With  Strong  Near-Surface  Structure  Complexity 


Network  Code 

Station  Code(s) 

G 

MBO 

GE 

APEZ,  CSS,  FODE,  KRIS,  KWP,  MELI, 
MRNI  NAI,  PUL,  RGN,  SFS,  SFUC,  SKD, 
VSU,  ZKR 

II 

ABKT,  GAR,  KIV 

IU 

TOL 

MN 

AIO,  AQU,  CEL,  CII,  CUC,  IDI,  MEB, 
TUE,  VAE 

SR 

MAIO 

XG 

BYKN,  DGSU,  EZRM,  HAMR,  HINS, 
HUB03  ,  KARS,  TTB01 

XI 

RIYD 

Using  &  Simplifying  the  Models 

To  some  extent  this  simplifying  assumption  allows  some  freedom  for  those  seeking 
to  use  the  models  as  a  starting  point.  Probably  the  best  constrained  parameters  in 
the  models  are  the  vertical  shear-wave  travel  times,  which  are  constrained  by  the  Ps 
wave  observed  on  most  of  the  receiver  functions.  Complexities  in  receiver  functions 
tend  to  introduce  more  roughness  into  the  structure,  in  an  effort  to  fit  small  but 
common  oscillations  following  the  P-wave.  These  arrivals  are  at  the  level  of  waves 
scattered  from  lateral  heterogeneity,  and  thus  inferences  drawn  on  one  dimensional 
structures  from  them  are  not  well  constrained.  They  remain  in  the  inversion  results 
because  the  simple  L2  norm  metric  of  misfit  does  not  distinguish  them  from  the 
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larger,  more  important  arrivals  that  are  reasonably  well  fit  (i.e.  the  residuals  of 
the  features  are  of  comparable  amplitudes).  Some  simplification  is  justified  when 
using  the  models  for  other  purposes.  How  much  a  model  can  be  adjusted  could 
be  tested  by  comparing  receiver  functions  and  dispersion  values  for  the  original  and 
modified  models.  Small  changes  in  the  models  may  not  strongly  affect  the  data  used 
to  constrain  the  model,  but  could  improve  predictions  of  other  data,  substantially. 
For  example,  even  a  minor  adjustment  of  the  wave  speed  gradient  beneath  the  crust 
could  strongly  affect  Pn  and  Sn  amplitudes.  Modifications  are  best  made  with  an 
eye  on  the  original  fits  to  the  data  and  with  some  appreciation  of  what  is  geologically 
reasonable  (Figure  3.6). 


Software  Availability 


An  implementation  of  the  methods  described  above  has  been  publicly  available  for 
about  two  years.  The  tools  are  integrated  into  Herrmann’s  Computer  Programs  in 
Seismology,  and  include  documentation. 
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Figure  3.6:  The  shear-speed  models  are  at  best  approximations  to  the  true  structure. 
Some  simplification  is  justified  when  using  the  models  for  other  purposes.  How 
much  a  model  can  be  adjusted  could  be  tested  by  comparing  receiver  functions  and 
dispersion  values  for  the  original  and  modified  models.  Modifications  are  best  made 
with  an  eye  on  the  original  fits  to  the  data  and  with  some  appreciation  of  what  is 
geologically  reasonable. 


Appendix  A 

Station  Summary  Tables 


In  this  appendix  we  include  summary  tables  that  list  the  locations,  number  of  ob¬ 
servations,  and  key  results  for  each  seismic  station  that  we  analyzed.  The  tables 
are  arranged  by  crustal  province,  and  each  station  is  listed  with  its  network  code 
first,  followed  by  the  station  code.  We  use  a  simple  crustal  classification  based  on 
the  map  found  in  Condie  [1993]  that  includes  shields,  platforms,  paleozoic  orogenic 
regions,  and  mesozoic-cenozoie  orogenic  regions.  Shields  are  platforms  are  generally 
similar,  except  for  a  sedimentary  veneer,  which  covers  the  basement  of  platforms. 
Each  result  is  ranked  in  the  range  from  A  to  D,  with  A  being  the  most  reliable 
results.  Specific  uncertainties  depend  on  the  particular  data  used  in  any  station 
analyses,  but  based  on  general  assumptions,  a  rule  of  thumb  uncertainty  for  crustal 
thickness  is  from  2  to  5  km.  We’ve  separated  out  three  stations  associated  with  the 
thick  crust  in  the  Pamir-Hindu  Kush  region.  For  information  on  Crust  2.0,  please 
visit  http://mahi.ucsd.edu/Gabi/rem.dir/crust/crust2.html. 

Table  A.l:  Shield  Station  Observations 


STATION 

LAT 

LON 

ELY 

(m) 

Data 

# 

RFTN 

# 

YEAR 

LOCATION 

G  TAM 
30<BAZ<80 

22.8 

5.53 

1377 

850 

132 

1990-04 

Africa 

G  TAM 

80<BAZ<360 

22.8 

5.53 

1377 

850 

132 

1990-04 

Africa 

GE  NAI 

-1.27 

36.8 

1992 

62 

23 

1995 

Africa 

GT  BGCA 

5.17 

18.42 

676 

353 

208 

1994-02 

Africa 

II  MBAR 

-6.01 

30.74 

1390 

208 

29 

1999-02 

Africa 

continued  on  next  page 
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Table  A.l:  continued 


STATION 

LAT 

LON 

ELV 

(m) 

Data 

# 

RFTN 

# 

YEAR 

LOCATION 

IU  FURI 

8.9 

36.69 

2545 

569 

37 

1997-02 

Africa 

IU  KMBO 

-1.12 

37.25 

1960 

452 

46 

1995-02 

Africa 

IU  KOWA 

14.5 

-4.01 

321 

39 

5 

1998-01 

Africa 

IU  MSKU 

-1.66 

13.61 

312 

99 

6 

1999-00 

Africa 

XD  AMBA 

-8.11 

33.26 

1415 

55 

8 

1994-95 

Africa-Tanzania 

XD  BASO 

-4.33 

35.14 

1694 

77 

12 

1994 

Africa- Tanzania 

XD  GOMA 

-4.84 

29.7 

880 

141 

25 

1994-95 

Africa-Tanzania 

XD  HALE 

-5.3 

38.62 

230 

198 

21 

1994-95 

Africa-Tanzania 

XD  INZA 

-5.2 

30.4 

975 

58 

6 

1994-95 

Africa-Tanzania 

XD  KIBA 

-5.32 

36.57 

1500 

155 

31 

1994-95 

A  fr  ica-Tanzania 

XD  KIBE 

-5.37 

36.47 

997 

228 

20 

1994-95 

Africa-Tanzania 

XD  KOMO 

-3.84 

36.72 

1114 

157 

24 

1994-95 

Africa-Tanzania 

XD  KOND 

-4.9 

35.8 

1419 

159 

14 

1994-95 

Africa-Tanzania 

XD  LONG 

-2.73 

36.7 

1380 

127 

22 

1994-95 

A  f  r  ica-  Tanz  an  i  a 

XD  MB  WE 

-4.96 

34.35 

1100 

136 

19 

1994-95 

A  frica-Tanzania 

XD  MITU 

-6.1 

34.06 

1566 

100 

18 

1994-95 

Africa-Tanzania 

XD  MTAN 

-7.91 

33.32 

1393 

105 

4 

1994 

Africa-Tanzania 

XD  MTOR 

-5.25 

35.4 

1100 

180 

33 

1994-95 

Africa-Tanzania 

XD  PAND 

-8.98 

33.24 

1248 

101 

12 

1994-95 

A  fr ica-Tanzan  ia 

XD  PUGE 

-4.71 

33.18 

1350 

113 

20 

1994-95 

Africa-Tanzania 

XD  RUNG 

-6.94 

33.52 

1230 

160 

25 

1994-95 

Africa-Tanzania 

XD  SING 

-4.64 

34.73 

1462 

187 

21 

1994-95 

Africa-Tanzania 

XD  TARA 

-3.89 

36.02 

1268 

135 

25 

1994-95 

Africa-Tanzania 

XD  TUND 

-9.3 

32.8 

1660 

126 

14 

1994-95 

Africa-Tanzania 

XD  URAM 

-5.08 

32.08 

1120 

127 

17 

1994-95 

Africa-Tanzania 

MN  WDD 

35.87 

41.52 

41 

474 

57 

1995-02 

Arabian-Plate 

IU  KONO 

59.65 

9.6 

216 

3473 

13 

1991-04 

Europe 

IU  KIEV 

50.7 

29.21 

163 

707 

194 

1995-03 

Europe 

XG  HUB03 

36.33 

74.66 

2469 

446 

48 

1999-01 

Middle-East 

XG  TRB01 

35.6 

75.22 

2072 

152 

5 

1996 

Middle-East 

XG  TTB01 

35.41 

74.59 

2610 

17 

5 

1995 

Middle-East 

XI  AFIF 

23.93 

43.04 

1116 

252 

32 

1995-97 

Middle-East 

XI  BISH 

19.92 

42.69 

1379 

27 

2 

1995 

Middle-East 

XI  HALM 

22.84 

44.32 

930 

343 

31 

1995-97 

Middle-East 

XI  RANI 

21.31 

42.78 

1001 

204 

18 

1995-97 

Middle-East 

XI  SODA 

18.29 

42.37 

2876 

289 

41 

1995-97 

Middle-East 

XI  TAIF 

21.28 

40.34 

2050 

43 

6 

1996 

Middle-East 

continued  on  next  page 
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Table  A.l:  continued 


STATION 

LAT 

LON 

ELV 

(m) 

Data 

# 

RFTN 

# 

YEAR 

LOCATION  | 

XI  UQSK 

25.78 

42.36 

950 

207 

11 

1996-97 

Middle-East 
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Table  A. 2:  Shield  Station  Results 


STATION 

H 

(km) 

Vp/Vs 

H(km) 
Crusl2. 0 

Vp/Vs 

Crust2.0 

TYPE 

Crust2.0 

RANK 

LOCATION 

G  TAM 
30<BAZ<8<3 

40 

1.6 

40 

1.77 

HI 

B-C 

Africa 

G  TAM 
80<BAZ<360 

33.8 

1.79 

40 

1.77 

HI 

B-C 

Africa 

GE  NAI 

39.2 

1.8 

38 

1.77 

X5 

B 

Africa 

GT  BGCA 

41.9 

1.71 

43 

1.77 

11 

A 

Africa 

II  MBAR 

34.7 

1.55 

39 

1.77 

G2 

C 

Africa 

IU  FURI 

38.6 

1.86 

31 

1.77 

G9 

C-D 

Africa 

IU  KMBO 

39.9 

1.76 

38 

1.77 

X5 

A 

Africa 

IU  KOWA 

42.6 

1.66 

41 

1.79 

G3 

C 

Africa 

IU  MSKU 

45.6 

1.78 

38 

1.79 

12 

C 

Africa 

XD  AMBA 

34.5 

1.76 

33 

1.77 

X2 

B 

Africa- Tanzania 

XD  BASO 

38.2 

1.67 

38 

1.77 

X5 

B-C 

Africa- Tanzania 

XD  GOMA 

43.3 

1.78 

43.5 

1.77 

GF 

A 

Africa- Tanzania 

XD  HALE 

35.2 

1.7 

41 

1.8 

13 

B 

Africa-Tanzania 

XD  INZA 

38.9 

1.77 

41.5 

1.76 

GE 

B 

Africa-Tanzania 

XD  KIBA 

39.2 

1.72 

36 

1.77 

IE 

A 

Africa-Tanzania 

XD  KIBE 

38.1 

1.7 

36 

1.77 

IE 

B 

Africa-Tanzania 

XD  KOMO 

37 

1.71 

38 

1.77 

X5 

A 

Africa-Tanzania 

XD  KOND 

38.3 

1.66 

38 

1.77 

X5 

A-B 

Africa-Tanzania 

XD  LONG 

38.5 

1.7 

38 

1.77 

X5 

A-B 

Africa-Tanzania 

XD  MBWE 

36.5 

1.75 

38 

1.77 

X5 

A 

Africa-Tanzania 

XD  MITU 

37.2 

1.74 

36 

1.77 

IE 

A 

Africa-Tanzania 

XD  MTAN 

37.9 

1.67 

41.5 

1.76 

GE 

B-C 

Africa-Tanzania 

XD  MTOR 

37.9 

1.74 

38 

1.77 

X5 

A 

Africa-Tanzania 

XD  PAND 

32.4 

1.8 

33 

1.77 

X2 

B 

Africa-Tanzania 

XD  PUGE 

35.7 

1.77 

37 

1.76 

Gl 

A 

Africa-Tanzania 

XD  RUNG 

39.2 

1.79 

41.5 

1.76 

GE 

B 

Africa-Tanzania 

XD  SING 

38.5 

1.78 

38 

1.77 

X5 

A-B 

Africa-Tanzania 

XD  TARA 

37.2 

1.75 

38 

1.77 

X5 

A-B 

Africa-Tanzania 

XD  TUND 

34.6 

1.7 

33 

1.77 

X2 

D 

Africa-Tanzania 

XD  URAM 

40.3 

2 

37 

1.76 

Gl 

C-D 

Africa-Tanzania 

MN  WDD 

35.7 

1.69 

36 

1.81 

ID 

B 

Arabian-Plate 

IU  KONO 

32.2 

1.71 

32 

1.82 

T5 

B-C 

Europe 

IU  KIEV 

39.5 

1.79 

41 

1.78 

DC 

B 

Europe 

XG  HUB03 

50 

1.66 

65 

1.76 

R2 

C-D 

Middle-East 

XG  TRB01 

71.6 

1.69 

65 

1.76 

R2 

C 

Middle-East 

XG  TTB01 

68.8 

1.68 

65 

1.76 

R2 

C 

Middle-East 

continued  on  next  page 


APPENDIX  A.  STATION  SUMMARY  TABLES 


48 


Table  A. 2:  continued 


STATION 

H 

(km) 

Vp/Vs 

H{km) 

Crust2.0 

Vp/Vs 
Crust2. 0 

TYPE 
Crust2. 0 

RANK 

LOCATION 

XI  AFIF 

36.4 

1.75 

43 

1.77 

11 

B 

Middle-East 

XI  BISH 

37.1 

1.82 

43 

1.77 

11 

B 

Middle-East 

XI  HALM 

34.6 

1.79 

41 

1.79 

D9 

B 

Middle-East 

XI  RANI 

38.1 

1.76 

43 

1.77 

11 

B 

Middle-East 

XI  SODA 

38.4 

1.74 

43 

1.77 

11 

B 

Middle-East 

XI  TAIF 

37.4 

1.74 

36 

1.77 

IE 

B 

Middle-East 

XI  UQSK 

36.5 

1.74 

43 

1.77 

11 

B 

Middle-East 
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Table  A. 3:  Platform  Stations  Observations 


STATION 

LAT 

LON~' 

ELV 

(m) 

Data 

# 

RFTN 

# 

YEAR 

LOCATION 

G  MBO 

14.39 

-16.96 

3 

81 

12 

2002-03 

Africa 

MN  GFA 

34.34 

9.1 

250 

134 

9 

1990 

1995-96 

Africa 

MN  KEG 

29.93 

31.83 

460 

651 

58 

1990-99 

Africa 

GE  SUW 

54.01 

23.18 

152 

1882 

4 

1995-04 

Europe 

G-GRC 

47.295 

3.074 

191 

336 

35 

2000-02 

Europe 

GE  TIRR 

44.46 

28.41 

77 

82 

11 

2003-04 

Europe 

GE  TRTE 

58.4 

26.7 

0 

654 

42 

1996-03 

Europe 

GE  VSU 

58.46 

26.73 

63 

223 

12 

2003-04 

Europe 

G  WUS 

41.1 

79.22 

1457 

2104 

423 

1990-02 

Middle-East 

GE  BGIO 

31.72 

35.09 

752 

246 

35 

1994-96 

Middle-East 

GE  EIL 

29.67 

34.95 

717 

210 

138 

2000-02 

Middle-East 

GE  JER 

31.77 

35.2 

770 

495 

89 

1996-03 

Middle-East 

GE  MHV 

54.96 

37.77 

150 

1669 

263 

1995-04 

Middle-East 

GE  PUL 

59.77 

30.32 

65 

969 

49 

1998-03 

Middle-East 

II  KIV 
0<BAZ<90 

43.96 

42.69 

1054 

1350 

284 

1990-02 

Middle- East 

II  KIV 
90<BAZ<360 

43.96 

42.69 

1054 

1350 

284 

1990-02 

Middle-East 

II  OBN 

55.11 

36.57 

160 

4512 

244 

1990-04 

Middle- East 

II  RAYN 

23.52 

45.5 

631 

581 

144 

1997-02 

Middle-East 

IU  GNI 

40.15 

44.74 

1609 

542 

68 

1991-02 

Middle-East 

KZ  AKT 

50.43 

58.02 

360 

98 

44 

1994-99 

Middle-East 

KZ  AKTK 

50.43 

58 

360 

702 

77 

1999-02 

Middle-East 

KZ  PDG 

43.32 

79.48 

1277 

331 

28 

1999-02 

Middle-East 

XI  RIYD 

24.722 

46.64 

717 

244 

13 

1996-97 

Middle-East 

XW  AHQI 

40.94 

78.45 

1959 

72 

29 

1999-00 

Middle-East 

APPENDIX  A.  STATION  SUMMARY  TABLES 


50 


Table  A. 4:  Platform  Station  Results 


STATION 

H 

(km) 

Vp/Vs 

H(km) 

Crust2.0 

Vp/Vs 

Crust.2.0 

TYPE 

Crust2.0 

RANK 

LOCATION 

G  MBO 

26.3 

1.98 

26 

1.84 

S8 

C 

Africa 

MN  GFA 

36.7 

1.89 

36 

1.8 

DR 

D 

Africa 

MN  KEG 

32.2 

1.64 

32 

1.81 

IH 

B-C 

Africa 

G-GRC 

29.3 

1.77 

31 

1.8 

M3 

B 

Europe 

GE  SUW 

41.2 

1.8 

46 

1.79 

DE 

C 

Europe 

GE  TIRR 

38.9 

1.68 

35 

1.8 

M4 

B 

Europe 

GE  TRTE 

44.4 

1.72 

41 

1.78 

DC 

B 

Europe 

GE  VSU 

42.6 

1.74 

41 

1.78 

DC 

C 

Europe 

G  WUS 

54.1 

1.71 

50 

1.76 

QE 

A-B 

Middle-East 

GE  BGIO 

33.9 

1.74 

32 

1.81 

T4 

B-C 

Middle- East 

GE  EIL 

29.5 

1.86 

32 

1.8 

T6 

B-C 

Middle-East 

GE  JER 

34.9 

1.72 

32 

1.81 

T4 

A-B 

Middle-East 

GE  MHV 

49.6 

1.61 

46.5 

1.79 

H4 

A-B 

Middle-East 

GE  PUL 

40.1 

1.96 

41 

1.78 

DC 

C-D 

Middle- East 

II  KIV 
0<BAZ<90 

38.3 

1.92 

43 

1.78 

PB 

C-D 

Middle-East 

II  KIV 

90<BAZ<360 

46 

1.74 

43 

1.78 

PB 

C-D 

Middle-East 

II  OBN 

50 

1.7 

46.5 

1.79 

H4 

B 

Middle-East- 

II  RAYN 

38 

1.77 

41 

1.79 

D9 

A 

Middle-East 

IU  GNI 

44.9 

1.71 

43 

1.78 

Q3 

B 

Middle- East 

KZ  AKT 

44.6 

1.97 

50 

1.76 

PE 

C-D 

Middle-East 

KZ  AKTK 

50.5 

1.86 

50 

1.76 

PE 

C 

Middle-East 

KZ  PDG 

44.4 

1.86 

50 

1.77 

P4 

B 

Middle-East 

XI  RIYD 

41.6 

1.64 

41 

1.79 

D6 

B 

Middle-East 

XW  AHQI 

57.3 

1.76 

50 

1.76 

PE 

B 

Mid  die- East 
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Table  A. 5:  Paleozoic  Stations  Observations 


STATION 

LAT 

LON 

ELV 

(m) 

Data 

# 

RFTN 

# 

YEAR 

LOCATION 

CZ  DPC 
0<BAZ<80 

50.36 

16.41 

760 

936 

170 

1992-03 

Europe 

CZ  DPC 
80<BAZ<360 

50.36 

16.41 

760 

936 

170 

1992-03 

Europe 

G  SSB 

45.28 

4.542 

700 

815 

285 

199002 

Europe 

GE  DSB 

53.25 

-6.38 

236 

1741 

89 

1993-04 

Europe 

GE  HLG 

54.19 

7.88 

41 

399 

31 

2001-04 

Europe 

GE  IBBN 

52.31 

7.76 

140 

840 

26 

1999-04 

Europe 

GE  KWP 

49.63 

22.71 

463 

1011 

50 

1999-04 

Europe 

GE  MORC 

49.78 

17.54 

740 

895 

275 

1993-02 

Europe 

GE  MTE 

40.4 

-7.54 

815 

402 

64 

1997-03 

Europe 

GE  RGN 

54.55 

13.32 

15 

1591 

134 

1995-04 

Europe 

GE  RUE 

52.48 

13.78 

40 

797 

78 

2000-04 

Europe 

GE  SELV 

37.24 

-3.73 

650 

257 

24 

2002-03 

Europe 

GE  SFS 

36.47 

-6.21 

21 

378 

9 

2001-04 

Europe 

GE  WLF 

49.7 

6.15 

295 

845 

211 

1994-03 

Europe 

GR  GRA1 

49.69 

11.22 

499.5 

2630 

4 

1990-04 

Europe 

II  ESK 

55.32 

-3.21 

242 

2364 

13 

1990-04 

Europe 

IU  GRFO 

49.7 

11.22 

425 

752 

194 

1994-02 

Europe 

IU  PAB 

39.54 

-4.35 

925 

825 

183 

1992-02 

Europe 

IU  TOL 

39.9 

-4.05 

480 

89 

11 

1991-92 

Europe 

II  AAK 
0<BAZ<80 

42.64 

74.49 

1645 

1200 

172 

1991-02 

Middle-East 

II  AAK 
80<BAZ<360 

42.64 

74.49 

1645 

1200 

172 

1991-02 

Middle-East 

II  GAR 

39 

70.32 

1300 

115 

23 

1990-92 

Middle- East 

KN  AML 

42.13 

73.69 

3400 

2298 

339 

1991-04 

Middle-East 

KZ  BRVK 

53.06 

70.28 

315 

2390 

514 

1994-04 

Middle-East 

KZ  CHK 

53.68 

70.6 

120 

690 

144 

1994-99 

Middle-East 

KZ  KKAR 

43.1 

70.5 

524.9 

217 

97 

2003-04 

Middle-East 

KZ  TLG 

43.23 

77.23 

1120 

1899 

162 

1994-02 

Middle-East 

KZ  VOS 

52.72 

70.98 

300 

1364 

330 

1994-04 

Middle-East 

KZ  ZRN 

52.95 

69 

380 

109 

38 

1994  &  1999 

Middle- East 
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Table  A. 6:  Paleozoic  Station  Results 


STATION 

H 

(km) 

Vp/Vs 

H(km) 

Crust2.0 

Vp/Vs 

Crust2.0 

TYPE 

Crust2.0 

RANK 

LOCATION 

CZ  DPC 
0<BAZ<80 

25.8 

1.86 

34 

1.8 

MD 

B-C 

Europe 

CZ  DPC 
80<BAZ<360 

29.8 

1.74 

34 

1.8 

MD 

B 

Europe 

G  SSB 

27.1 

1.69 

31 

1.8 

Ml 

B 

Europe 

GE  DSB 

30.5 

1.7 

31 

1.79 

M5 

B 

Europe 

GE  HLG 

28.8 

1.74 

31 

1.82 

C4 

C 

Europe 

GE  IBBN 

30.9 

1.76 

31 

1.8 

Ml 

B-C 

Europe 

GE  KWP 

35.8 

1.73 

46 

1.79 

DE 

B 

Europe 

GE  MORC 

31.6 

1.79 

31 

1.79 

M5 

B 

Europe 

GE  MTE 

29 

1.76 

31 

1.79 

M5 

B 

Europe 

GE  RGN 

29.6 

1.73 

31 

1.8 

Ml 

B 

Europe 

GE  RUE 

29.2 

1.73 

31 

1.8 

MH 

B 

Europe 

GE  SELV 

41 

1.84 

32 

1.81 

T4 

C-D 

Europe 

GE  SFS 

32 

1.93 

32 

1.81 

T4 

D 

Europe 

GE  WLF 

28.8 

1.78 

31 

1.8 

M3 

B 

Europe 

GR  GRA1 

29.2 

1.83 

31 

1.8 

M3 

C-D 

Europe 

II  ESK 

28.2 

1.82 

31 

1.79 

M5 

C 

Europe 

IU  GRFO 

29.7 

1.76 

31 

1.8 

M3 

A 

Europe 

IU  PAB 

31.5 

1.7 

31 

1.79 

M5 

A 

Europe 

IU  TOL 

33.2 

1.73 

31 

1.79 

M5 

B 

Europe 

II  AAK 
0<BAZ<80 

50 

J 

1.78 

50 

1.77 

P4 

B 

Middle-East 

II  AAK 

80<BAZ<360 

53.6 

1.63 

50 

1.77 

P4 

B 

Middle- East 

II  GAR 

54.5 

1.75 

50 

1.77 

P4 

C 

Middle- East 

KN  AML 

53.3 

1.93 

45 

1.77 

PC 

C-D 

Middle-East 

KZ  BRVK 

47.2 

1.73 

41 

1.78 

DC 

B 

Middle-East 

KZ  CHK 

41.2 

1.72 

41 

1.78 

DC 

B 

Middle-East 

KZ  KKAR 

40.2 

1.65 

45 

1.77 

PC 

B 

Middle-East 

KZ  TLG 

46.4 

1.72 

50 

1.77 

P4 

C 

Middle-East 

KZ  VOS 

38.2 

2.1 

41 

1.78 

DC 

C-D 

Middle-East 

KZ  ZRN 

44.5 

1.74 

41 

1.79 

D9 

B 

Middle-East 
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Table  A. 7:  Mesozoic-Cenozoic  Stations  Observations 


STATION 

LAT 

LON 

ELV 

(m) 

Data 

# 

RFTN 

# 

YEAR 

LOCATION 

GE  MELI 

35.29 

-2.94 

40 

732 

18 

1999-04 

Africa 

MN  MDT 

32.82 

-4.61 

1200 

159 

34 

1990-91 

1994-95 

Africa 

MN  MEB 

36.3 

2.73 

500 

106 

13 

1992-93 

Africa 

G  ATD 

11.53 

42.85 

610 

780 

192 

1993-02 

Africa-Rift 

G  ECH 

48.22 

7.15 

580 

874 

194 

1990-02 

Europe 

GE  APE 

37.07 

25.53 

620 

596 

41 

2000-04 

Europe 

GE  APEZ 

34.98 

24.89 

435 

528 

30 

2000-03 

Europe 

GE  CART 

37.59 

-1 

65 

1042 

56 

1997-04 

Europe 

GE  FODE 

35.38 

24.96 

50 

556 

41 

2000-03 

Europe 

GE  GVD 

34.84 

24.09 

180 

731 

7 

1999-03 

Europe 

GE  KRIS 

35.18 

25.5 

850 

898 

55 

1998-04 

Europe 

GE  MAHO 

39.9 

4.27 

15 

738 

32 

1999-04 

Europe 

GE  MLR 

45.5 

25.95 

1378 

678 

53 

1994-02 

Europe 

GE  PSZ 

47.92 

19.9 

549 

1046 

178 

1995-04 

Europe 

GE  SANT 

36.37 

25.46 

540 

942 

11 

1998-04 

Europe 

GE  SFUC 

36.64 

-6.17 

88 

406 

38 

199601 

Europe 

GE  SKD 

35.41 

23.93 

306 

950 

83 

1998-04 

Europe 

GE  STU 

48.77 

9.2 

360 

842 

218 

1994-03 

Europe 

GE  ZKR 

35.11 

26.22 

270 

87 

4 

2003-04 

Europe 

II  BFO 

48.33 

8.33 

589 

551 

83 

1996-02 

Europe 

MN  AIO 

37.97 

15.23 

751.4 

248 

2 

2001-03 

Europe 

MN  AQU 

42.35 

13.41 

710 

1033 

140 

19903 

Europe 

MN  BGY 

44.8 

20.52 

250 

224 

28 

1991-95 

1997-98 

Europe 

MN  BNI 

45.1 

6.68 

1395 

549 

48 

1990-92 

1996-03 

Europe 

MN  CEL 

38.26 

15.89 

702 

50 

5 

2003-04 

Europe 

MN  CII 

41.72 

14.03 

910 

405 

7 

1994-96 

1999-02 

Europe 

MN  CLTB 

37.58 

13.22 

949 

274 

16 

2001-04 

Europe 

MN  CUC 

39.99 

15.81 

637 

81 

9 

2003-04 

Europe 

MN  IDI 

35.29 

24.89 

750 

524 

37 

1996-04 

Europe 

MN  TRI 

45.71 

13.25 

161 

222 

71 

1996-99 
&  2003 

Europe 

MN  TUE 

46.47 

9.35 

1924 

272 

2 

2001-04 

Europe 

MN  VAE 

37.47 

14.35 

735.1 

280 

1 

2001-03 

Europe 

MN  VSL 

39.5 

9.4 

370 

1001 

229 

1990-03 

Europe 
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Table  A. 7:  continued 


STATION 

LAT 

LON 

ELV 

(m) 

Data 

# 

RFTN 

# 

YEAR 

LOCATION 

MN  VTS 

42.62 

2.24 

1490 

391 

141 

1996-02 

Europe 

AS  KAAO 

34.54 

69.04 

19 

1079 

72 

1977-82 

Middle-East 

GE  CSS 

34.96 

33.33 

396 

1026 

221 

1999-04 

Middle-East 

GE  ISP 

37.84 

30.51 

1100 

500 

113 

1996-03 

Middle-East 

GE  KSDI 

33.19 

35.7 

123 

402 

31 

2002-04 

Middle-East 

GE  MALT 

38.31 

38.43 

1120 

714 

160 

2000-04 

Middle-East 

II  ABKT 

37.93 

58.12 

678 

622 

293 

1993-02 

Middle- East 

II  NIL 

33.65 

73.27 

629 

652 

23 

1994-00 

Middle-East 

IL  IR1 

35.42 

50.69 

200 

209 

50 

1976-79 

Middle-East 

IL  IR2 

35.67 

50.9 

200 

209 

18 

1976-79 

Middle- East 

IL  IR3 

35.48 

51.02 

200 

209 

40 

1976-79 

Middle-East 

IL  IR4 

35.24 

50.9 

200 

209 

17 

1976-79 

Middle-East 

IL  IR5 

35.21 

50.58 

200 

209 

12 

1976-79 

Middle- East 

IL  IR6 

35.47 

50.43 

200 

209 

7 

1976-79 

Mid  die- East 

IL  IR7 

35.7 

50.61 

200 

209 

33 

1976-79 

Middle-East 

IU  ANTO 

39.87 

32.79 

883 

376 

70 

1992-98 

Mid  die- East 

SR  MAIO 

36.3 

59.5 

1000 

572 

84 

1975-78 

Middle-East 

XG  AGIN 

38.94 

38.71 

945 

356 

30 

1999-01 

Mid  die- East 

XG  AHLT 

38.75 

42.48 

1738 

282 

34 

1999-01 

Middle- East 

XG  BNGL 

38.9 

40.6 

1160 

440 

72 

1999-01 

Middle-East 

XG  BTLS 

38.43 

42.12 

1722 

327 

43 

1999-01 

Middle-East 

XG  BYBT 

40.23 

40.26 

1684 

301 

29 

1999-01 

Middle- East 

XG  BYKN 

38.17 

41.78 

816 

230 

41 

1999-01 

Middle-East 

XG  CMCY 

39.92 

43.2 

1995 

295 

34 

1999-01 

Middle-East 

XG  DGRL 

41.06 

43.33 

2027 

358 

48 

1999-01 

Middle-East 

XG  DGSU 

39.13 

42.73 

1646 

175 

18 

1999-01 

Middle-East 

XG  DYBR 

37.82 

40.32 

619 

337 

38 

1999-01 

Middle-East 

XG  ERGN 

38.26 

39.73 

922 

445 

47 

1999-01 

Middle-East 

XG  EZRM 

40.1 

41.36 

1910 

442 

29 

1999-01 

Middle- East 

XG  HAMR 

39.61 

42.99 

1714 

239 

26 

1999-01 

Middle-East 

XG  HINS 

39.35 

41.7 

1727 

317 

70 

1999-01 

Middle- East 

XG  HRPT 

38.7 

39.25 

142 

429 

78 

1999-01 

Middle-East 

XG  HRSN 

39.94 

42.29 

1800 

446 

32 

1999-01 

Middle-East 

XG  ILIC 

39.45 

38.57 

1300 

223 

50 

1999-01 

Middle-East 

XG  IMRL 

39.88 

38.12 

1708 

172 

7 

1999-01 

Middle-East 

XG  KARS 

40.62 

43.07 

1904 

190 

39 

1999-01 

Middle-East 

XG  KOTK 

40.22 

43 

1396 

247 

55 

1999-01 

Middle-East. 
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Table  A. 7:  continued 


STATION 

LAT 

LON 

ELV 

(m) 

Data 

# 

RFTN 

# 

YEAR 

LOCATION 

XG  KRLV 

39.37 

40.99 

1894 

381 

34 

1999-01 

Middle-East 

XG  KTLN 

37.95 

41.7 

791 

228 

20 

1999-01 

Middle-East 

XG  KYPR 

37.56 

41.17 

1179 

274 

62 

1999-01 

Middle-East 

XG  MRDN 

37.29 

40.7 

674 

425 

61 

1999-01 

Middle-East 

XG  MSDY 

40.46 

37.78 

1152 

458 

52 

1999-01 

Middle-East 

XG  MUSH 

38.76 

41.48 

1554 

317 

49 

1999-01 

Middle-East 

XG  SILN 

38.14 

41.04 

917 

237 

38 

1999-01 

Mid  die- East 

XG  SIRN 

40.2 

39.12 

1496 

445 

62 

1999-01 

Middle- East 

XG  UZML 

39.71 

39.72 

1593 

579 

2 

1999-01 

Middle-East 
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Table  A. 8:  Mesozoic-Cenozoic  Station  Results 


STATION 

H 

(km) 

Vp/Vs 

H(km) 

Crust2.0 

Vp/Vs 
Crust2. 0 

TYPE 

Crust2.0 

RANK 

LOCATION 

GE  MELI 

34.7 

1.52 

32 

1.81 

T4 

C 

Africa 

MN  MDT 

33.2 

1.72 

36 

1.77 

IE 

B 

Africa 

MN  MEB 

26.2 

1.9 

24.5 

1.86 

S3 

C 

Africa 

G  ATI) 

22.6 

1.73 

17 

1.79 

X4 

B-C 

Africa-Rift 

G  ECH 

25.1 

1.72 

31 

1.8 

M3 

C 

Europe 

GE  APE 

21.5 

1.75 

26.5 

1.86 

S4 

B 

Europe 

GE  APEZ 

29.8 

2.1 

27 

1.86 

S5 

D 

Europe 

GE  CART 

20.6 

1.79 

24.5 

1.86 

S6 

C 

Europe 

GE  FODE 

31.7 

1.7 

27 

1.86 

S5 

C 

Europe 

GE  GVD 

26 

1.77 

27 

1.86 

S5 

C 

Europe 

GE  KRIS 

28.2 

1.81 

27 

1.86 

S5 

C 

Europe 

GE  MAHO 

19.6 

1.7 

18 

1.89 

B8 

B 

Europe 

GE  MLR 

34.26 

1.64 

34.5 

1.8 

MD 

C 

Europe 

GE  PSZ 

27.5 

1.75 

27.5 

1.8 

ME 

B 

Europe 

GE  SANT 

26.5 

1.81 

26.5 

1.86 

S4 

C 

Europe 

GE  SFUC 

32.25 

1.64 

32 

1.81 

T4 

C 

Europe 

GE  SKD 

37.5 

1.96 

26 

1.84 

S8 

D 

Europe 

GE  STU 

24.3 

1.75 

27.5 

1.79 

MG 

A 

Europe 

GE  KRIS 

28.2 

1.81 

27 

1.86 

S5 

C 

Europe 

GE  ZKR 

32.9 

2 

27 

1.86 

S5 

C 

Europe 

II  BFO 

26.1 

1.66 

27.5 

1.8 

MG 

B 

Europe 

MN  AIO 

31.2 

1.99 

26.5 

1.86 

S4 

D 

Europe 

MN  AQU 

37.52 

1.7 

31 

1.79 

M5 

B 

Europe 

MN  BGY 

27.12 

1.78 

31 

1.79 

M5 

B 

Europe 

MN  BNI 

44.9 

1.86 

38 

1.77 

PA 

B 

Europe 

MN  CEL 

27 

1.83 

24.5 

1.86 

S6 

B-C 

Europe 

MN  CII 

32 

1.84 

31 

1.8 

Ml 

C-D 

Europe 

MN  CLTB 

30.1 

1.78 

26 

1.83 

S7 

C 

Europe 

MN  CUC 

17.9 

1.91 

24.5 

1.86 

S6 

C 

Europe 

MN  IDI 

28.8 

1.68 

27 

1.86 

S5 

C 

Europe 

MN  TRI 

36 

1.77 

37 

1.8 

M8 

B 

Europe 

MN  TUE 

30.1 

1.88 

38 

1.77 

PA 

C 

Europe 

MN  VAE 

26.5 

1.79 

26.5 

1.86 

S4 

C-D 

Europe 

MN  VSL 

30.6 

1.69 

22 

1.82 

SI 

B 

Europe 

MN  VTS 

38.1 

1.73 

36 

1.79 

MG 

B 

Europe 

AS  KAAO 

45 

1.89 

38 

1.79 

DT 

D 

Middle-East 

GE  CSS 

35.9 

1.91 

24.5 

1.86 

S6 

C-D 

Middle-East 
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Table  A. 8:  continued 


STATION 

H 

(km) 

Vp/Vs 

H(km.) 
Crust  2.0 

Vp/Vs 
Crust2. 0 

TYPE 

Crust2.0 

RANK 

LOCATION 

GE  ISP 

36 

1.68 

38 

1.78 

P2 

B-C 

Middle-East 

GE  KSDI 

30.5 

1.89 

26.5 

1.86 

S4 

C 

Middle-East 

GE  MALT 

37.3 

1.54 

43 

1.77 

P6 

B 

Middle-East 

II  ABKT 

41.8 

1.7 

46 

1.79 

DE 

B-C 

Middle- East 

II  NIL 

51 

1.76 

50 

1.78 

PP 

B 

Middle- East 

IL  IR1 

48 

1.77 

41 

1.79 

D9 

C 

Middle-East 

IL  IR2 

43.2 

1.82 

41 

1.79 

D9 

C 

Mid  die- East. 

IL  IR3 

45.18 

1.8 

41 

1.79 

D9 

C 

Middle-East 

IL  IR4 

42 

1.96 

41 

1.79 

D9 

D 

Middle- East 

IL  IR5 

49 

1.85 

41 

1.79 

D9 

D 

Middle-East 

IL  IR6 

44.7 

1.87 

41 

1.79 

D9 

D 

Middle-East 

IL  IR7 

41.5 

1.88 

41 

1.79 

D9 

C 

Middle- East 

IU  ANTO 

35 

1.71 

46 

1.77 

PI 

B-C 

Middle-East 

SR  MAIO 

48.43 

1.77 

46 

1.79 

DE 

C 

Middle- East 

XG  AGIN 

37.3 

1.76 

43 

1.77 

P6 

B 

Middle- East 

XG  AHLT 

42.9 

1.76 

43 

1.77 

P6 

B-C 

Middle- East 

XG  BNGL 

42.8 

1.84 

43 

1.77 

P6 

B-C 

Middle-East 

XG  BTLS 

40.5 

1.83 

43 

1.77 

P6 

C 

Middle-East 

XG  BYBT 

38.8 

1.89 

36 

1.  79 

DQ 

C 

Middle- East 

XG  BYKN 

46.2 

1.63 

43 

1.77 

P6 

B-C 

Middle-East 

XG  CMCY 

41.7 

1.83 

43 

1.77 

P6 

D 

Middle- East 

XG  DGRL 

37 

1.86 

39 

1.76 

Q2 

C-D 

Middle-East 

XG  DGSU 

41 

1.84 

43 

1.77 

P6 

C 

Middle-East 

XG  DYBR 

35.3 

1.67 

36 

1.78 

Q5 

C 

Middle-East 

XG  ERGN 

45.2 

1.7 

43 

1.77 

P6 

C-D 

Middle-East 

XG  EZRM 

40.3 

1.79 

36 

1.79 

DQ 

C-D 

Middle- East 

XG  HAMR 

45.6 

1.69 

43 

1.77 

P6 

C 

Middle-East 

XG  HINS 

44.2 

1.57 

43 

1.77 

P6 

C 

Mid  die- East 

XG  HRPT 

35.6 

1.81 

43 

1.77 

P6 

CD 

Mid  die- East 

XG  HRSN 

41 

1.81 

43 

1.77 

P6 

C 

Middle-East 

XG  ILIC 

42 

1.91 

43 

1.77 

P6 

C 

Middle-East 

XG  IMRL 

43.3 

1.84 

43 

1.77 

P6 

C 

Middle-East 

XG  KARS 

36.5 

2.17 

39 

1.76 

Q2 

D 

Middle-East 

XG  KOTK 

40.9 

1.83 

39 

1.76 

Q2 

C 

Middle-East 

XG  KRLV 

40.4 

1.83 

43 

1.77 

P6 

C 

Middle-East 

XG  KTLN 

24.7 

2.19 

36 

1.78 

Q5 

D 

Middle-East 

XG  KYPR 

32.6 

1.85 

36 

1.78 

Q5 

C 

Middle-East 
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Table  A. 8:  continued 


STATION 

H 

(km) 

Vp/Vs 

H{km) 

Crust2.0 

Vp/Vs 
Crust2. 0 

TYPE 
Crust2. 0 

RANK 

LOCATION 

XG  MRDN 

35.9 

1.76 

36 

1.78 

Q5 

B 

Middle-East, 

XG  MSDY 

38.1 

1.84 

41 

1.79 

D6 

C-D 

Middle-East 

XG  MUSH 

44.7 

1.75 

43 

1.77 

P6 

C 

Mid  die- East 

XG  SILN 

41.3 

1.64 

43 

1.77 

P6 

C 

Middle-East 

XG  SIRN 

38.5 

1.68 

38 

1.8 

D3 

C 

Middle- East 

XG  UZML 

40.4 

1.79 

43 

1.77 

P6 

D 

Middle- East 
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Table  A. 9:  Pamir  Stations  Observations 


STATION 

LAT 

LON 

ELV 

(m) 

Data 

# 

RFTN 

# 

YEAR 

LOCATION 

XG  HUB03 

36.33 

74.66 

2469 

446 

75 

1999-01 

Middle-East 

XG  TRB01 

35.6 

75.22 

2072 

152 

5 

1996 

Middle-East 

XG  TTB01 

35.41 

74.59 

2610 

17 

5 

1995 

Middle-East 

Table  A. 10:  Pamir  Station  Results 


STATION 

H 

(km) 

Vp/Vs 

H(km ) 
Crust2.0 

Vp/Vs 

Crust2.0 

TYPE 
Crust  2.0 

RANK 

LOCATION 

XG  HUB03 

50 

1.66 

65 

1.76 

R2 

C-D 

Middle-East 

XG  TRB01 

71.6 

1.69 

65 

1.76 

R2 

C 

Middle-East 

XG  TTB01 

68.8 

1.68 

65 

1.76 

R2 

C 

Middle-East 

Appendix  B 

Inversion  Summaries 


In  this  appendix  we  include  summaries  for  each  of  the  171  stations  for  which  we  have 
attempted  a  shear-wave  speed  inversion.  Please  note,  this  is  a  complete  summary 
and  we  do  not  believe  that  all  observations  could  be  explained  with  these  relatively 
simple  isotropic  plane-layered  earth  models.  These  data  generally  result  in  wild 
shear-wave  speed  models,  that  are  unacceptable  in  terms  of  simplicity  or  any  under¬ 
standing  of  rock  physics.  Data  of  quality  D  in  the  summary  tables  in  the  previous 
appendix  identify  stations  that  we  think  were  too  complicated  or  had  too  little  data 
to  produce  a  reliable  result.  Each  page  summarizes  the  results  at  one  station.  The 
upper  left  shows  some  summary  results  from  the  H-K  stacking,  and  the  azimtuh 
and  ray  parameter  ranges.  Receiver  functions  from  two  band  widths  are  included  in 
each  inversion.  The  red  curves  correspond  an  analysis  based  on  the  Harvard  group 
velocity  model,  the  blue  to  the  LLNL  group  velocity  model.  The  gray  region  shows 
the  range  of  typical  crustal  shear-wave  speeds  reported  in  Christensen  and  Mooney 
[1995]. 
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Station:  AS  KAAO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.059  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  51,  G2.50  =  56 
Crustal  Thickness 
HK  method  =  45  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.89 
Crust2 . 0  =  1.785 


-5  0  5  10  15  20  25  30  35  40  45  50  55  60 

Time  (s) 


VS  (km/s) 


Period  (s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 
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Station:  CZ  DPC 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0517  (s/km) 

Influence  factor  =  0.5 

Smoothing  factor  =  0.5 

Number  of  RFTN 

G1.00  =  195.5,  G2.50  =  160 

Crustal  Thickness 

HK  method  =  28  km 

Crust2.0  =  34  km 

Vp/Vs  Ratio 

HK  method  =  0 

Crust2 . 0  =  1.8 


Time  (s) 


VS  (km/s) 


0  20  40  60  80  100  120  140  160  180  200 

Period  (s) 
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Station:  GE  APE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0618  (s/km) 
Influence  factor  -  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  41,  G2.50  =  12 
Crustal  Thickness 
HK  method  =  21.5  km 
Crust2.0  =  26.5  km 
Vp/Vs  Ratio 
HK  method  =  1.75 
Crust2 . 0  =  1.86 
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Station:  GE  APEZ 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0547  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  13,  G2.50  =  3 
Crustal  Thickness 
HK  method  =  29.8  km 
Crust2.0  =  27  km 
Vp/Vs  Ratio 
HK  method  =2.1 
Crust2.0  =  1.86 
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Group  Velocity  (km/s)  Receiver  Function  (1/s) 
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Station:  GE  BGIO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0506  (s/kro) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  51,  G2.50  =  38 
Crustal  Thickness 
HK  method  =  33.9  km 
Crust2.0  =  32  km 
Vp/Vs  Ratio 
HK  method  =  1.74 
Crust2 .0  =  1.81 
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Station:  GE  CSS 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.049  ( s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  100,  G2.50  =  40 
Crustal  Thickness 
HK  method  =  36  km 
Crust2 . 0  =  24.5  km 
Vp/Vs  Ratio 
HK  method  =  1.91 
Crust2 .0  =  1.86 
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Station:  GE  DSB 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0511  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  66,  G2.50  =  46 
Crustal  Thickness 
HK  method  =  30.5  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.7 
Crust2.0  =  1.79 
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Station:  6E  EIL 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0524  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00' =  37,  G2.50  =  57 
Crustal  Thickness 
HK  method  =  29.5  km 
Crust2.0  =  32  km 
Vp/Vs  Ratio 
HK  method  =  1.86 
Crust2 .0  =  1.8 
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Station:  GE  FODE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0466  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  12,  G2.50  =  7 
Crustal  Thickness 
HK  method  =  31.7  km 
Crust2 .0  =  27  km 
Vp/Vs  Ratio 
HK  method  =  1.7 
Crust2 .0  -  1.86 
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Station:  GE  GVD 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0515  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  4,  G2.50  =  4 
Crustal  Thickness 
HK  method  =  26  km 
Crust2.0  =  27  km 
Vp/Vs  Ratio 
HK  method  -  1.77 
Crust2 .0  -  1.86 
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Group  Velocity  (km/s)  Receiver  Function  (1/s) 
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Station:  6E  HLG 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0589  (s/km) 
Influence  factor  =0*5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  17 #  G2.50  =  15 
Crustal  Thickness 
HK  method  =  28.8  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.74 
Crust2 . 0  =  1.82 
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Period  (s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 
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Station:  GE  IBBN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.057  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  14,  G2.50  =  16 
Crustal  Thickness 
HK  method  =  30.9  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 . 0  =  1.86 
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Station:  GE  ISP 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.052  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  55 ,  G2.50  =  55 
Crustal  Thickness 
HK  method  =  36  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.68 
Crust2.0  =1.78 
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Station:  GE  JER 


Shear  Velocity  Models 


Back  a2imuth  range  =  All 
Average  ray  parameter 
P  =  0.053  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  48,  G2.50  =  41 
Crustal  Thickness 
HK  method  =  34.3  km 
Crust2.0  =  32  km 
Vp/Vs  Ratio 
HK  method  =  1.74 
Crust2 . 0  =  1.81 
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Station:  GE  KRIS 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0521  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  33,  G2.50  =  27 
Crustal  Thickness 
HK  method  =  28.2  km 
Crust2.0  =  27  km 
Vp/Vs  Ratio 
HK  method  =  1.81 
Crust2 .0  =  1.86 
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Station:  GE  KSDI 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0529  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  »  18,  G2.50  =  15 
Crustal  Thickness 
HK  method  =  30.45  km 
Crust2 . 0  =  26.5  km 
Vp/Vs  Ratio 
HK  method  =  1.89 
Crust2 .0  =  1.86 
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Group  Velocity  (km/s)  Receiver  Function  (1/s) 
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Station:  GE  KWP 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0498  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  29,  G2.50  =  26 
Crustal  Thickness 
HK  method  =  35.8  km 
Crust2.0  =  46  km 
Vp/Vs  Ratio 
HK  method  =  1.73 
Crust2 . 0  =  1.79 
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Station:  GE  MAHO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.051  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  40,  G2.50  =  26 
Crustal  Thickness 
HK  method  =  19.6  km 
Crust2.0  =  18  km 
Vp/Vs  Ratio 
HK  method  =  1.70 
Crust2 . 0  =  1.89 
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Station:  GE_MALT 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0534  (s/km) 

Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  74,  G2.50  =  97 
Crustal  Thickness 
HK  method  =  37.29  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.54 
Crust2 . 0  =1.77 


Shear  Velocity  Models 
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Station:  GE_MELI 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.058  (s/km) 

Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  8 ,  G2.50  =  5 
Crustal  Thickness 
HK  method  =  34.7  km 
Crust2.0  =  32  km 
Vp/Vs  Ratio 
HK  method  =  1.52 
Crust2.0  =  1.81 


Shear  Velocity  Models 
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Station:  GE  MHV 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  *  0.057  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  179,  G2.50  =  158 
Crustal  Thickness 
HK  method  -  49.6  km 
Crust2.0  =  46.5  km 
Vp/Vs  Ratio 
HK  method  =  1.61 
Crust2 . 0  =  1.79 
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Station:  6E  MLR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0561  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  46,  G2.50  =  41 
Crustal  Thickness 
HK  method  =  34.3  km 
Crust2 .0  =  34.5  km 
Vp/Vs  Ratio 
HK  method  =  1.64 
Crust2 . 0  -  1.8 
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Station:  GE  MORC 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0517  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  227,  G2.50  =  225 
Crustal  Thickness 
HK  method  =  31.6  km 
Crust2.0  =  3 1  km 
Vp/Vs  Ratio 
HK  method  =  1.79 
Crust2 . 0  =  1.79 
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Station:  GE_MRNI 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0534  (s/km) 

Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  50,  G2.50  =  42 
Crustal  Thickness 
HK  method  =  33.1  km 
Crust2.0  =  26.5  km 
Vp/Vs  Ratio 
HK  method  =  1.63 
Crust2 . 0  =  1.86 
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Station:  GE  MTE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0535  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  51.5,  G2.50  =  41 
Crustal  Thickness 
HK  method  =  29  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 . 0  =  1.79 


-5  0  5  10  15  20  25  30  35  40  45  50  55  60 

Time  (s) 


VS  (km/s) 


0  20  40  60  80  100  120  140  160  180  200 

Period  (s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 


APPENDIX  B.  INVERSION  SUMMARIES 


86 


Station:  GE  NAI 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0484  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  8,  G2.50  =  9 
Crustal  Thickness 
HK  method  =  39.2  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.8 
Crust2 .0  =  1.77 
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Station:  GE  PSZ 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0572  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  60 ,  G2.50  =  101 
Crustal  Thickness 
HK  method  =  27.5  km 
Crust2 . 0  =  27.5  km 
Vp/Vs  Ratio 
HK  method  =  1.75 
Crust2 .0  =  1.8 
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Station:  GE  PUL 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.06  (s/km) 

Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  43,  G2.50  =  42 
Crustal  Thickness 
HK  method  =  40.1  km 
Crust2.0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.96 
Crust2.0  =  1.78 


1 1 1 1  ’ 1 1 m  ■  i i ' "  >  . ■  i 1 1 1 1 1 1  1  ■  i ' ' "  i  1 1  ’ r r 

-5  0  5  10  15  20  25  30  35  40  45  50  55  60 

Time  (s) 


VS  (km/s) 


Period  (s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 


APPENDIX  B.  INVERSION  SUMMARIES 


89 


Station:  6E  RGN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0559  (s/km) 
Influence  factor  -  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  65,  G2.50  =  34 
Crustal  Thickness 
HK  method  =  29.7  km 
Crust2 . 0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.73 
Crust2 . 0  =  1.8 
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Station:  GE  RUE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0534  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  20,  G2.50  =  30 
Crustal  Thickness 
HK  method  =  29.2  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.73 
Crust2 . 0  =1.8 
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Station:  GE  SANT 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0613  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  5,  G2.50  =  3 
Crustal  Thickness 
HK  method  =  26.5  km 
Crust2.0  =  26.5  km 
Vp/Vs  Ratio 
HK  method  =  1.81 
Crust2 . 0  =  1.86 
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Station:  GE  SELV 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0523  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  9,  G2.50  =  6 
Crustal  Thickness 
HK  method  =  41  km 
Crust2 .0  =  32  km 
Vp/Vs  Ratio 
HK  method  =  1.84 
Crust2.0  =  1.81 
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Station:  GE  SFS 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.061  { s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  3,  G2.50  =  5 
Crustal  Thickness 
HK  method  =  31.6  km 
Crust2.0  =  32  km 
Vp/Vs  Ratio 
HK  method  =  1.95 
Crust2 .0  =  1.81 
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Station:  GE_SFUC 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0498  (s/km) 

Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  19,  G2.50  =  12 
Crustal  Thickness 
HK  method  =  32.2  km 
Crust2.0  =  32  km 
Vp/Vs  Ratio 
HK  method  =  1.64 
Crust2.0  =  1.81 
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Station:  6E  SKD 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0559  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  18,  G2.50  -  14 
Crustal  Thickness 
HK  method  =  37.5  km 
Crust2.0  =  26  km 
Vp/Vs  Ratio 
HK  method  =  1.96 
Crust2.0  =  1.84 
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Station:  GE  STU 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0519  (s/km) 

Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  163,  G2.50  =  164 
Crustal  Thickness 
HK  method  =  24.3  km 
Crust2.0  =  27.5  km 
Vp/Vs  Ratio 
HK  method  =  1.75 

Crust2 . 0  =  1.79  ^ 
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Station:  GE  SUW 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.053  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  1,  G2.50  ==  1 
Crustal  Thickness 
HK  method  =  41.2  km 
Crust2 .0  =  46  km 
Vp/Vs  Ratio 
HK  method  =  1.8 
Crust2 . 0  =  1.79 
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Station:  GE  TIRR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0477  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  7,  G2.50  =  3 
Crustal  Thickness 
HK  method  =  38.9  km 
Crust2.0  =  35  km 
Vp/Vs  Ratio 
HK  method  =  1.68 
Crust2 . 0  =1.8 
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Station:  GE  TRTE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
p  =  0.0583  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  32,  G2.50  =  39 
Crustal  Thickness 
HK  method  =  44.4  km 
Crust2.0  =  4 1  km 
Vp/Vs  Ratio 
HK  method  =  1.72 
Crust2.0  =  1.78 
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Station:  GE  VSU 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0627  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  11,  G2.50  =  11 
Crustal  Thickness 
HK  method  =  42.6  km 
Crust2.0  -  41  km 
Vp/Vs  Ratio 
HK  method  =  1.741 
Crust2 . 0  =  1.78 
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Station:  GE  WLF 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.051  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  175,  G2.50  =  155 
Crustal  Thickness 
HK  method  =  28.8  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.77 
Crust2 .0  =  1.8 
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Station:  GE  ZKR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.042  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  4 ,  G2.50  =  1 
Crustal  Thickness 
HK  method  =  32.9  km 
Crust2.0  =  27  km 
Vp/Vs  Ratio 
HK  method  =  2 
Crust2 .0  =  1.86 
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Station:  GR  GRAl 


Shear  Velocity  Models 


Back  azimuth  range  -  All 
Average  ray  parameter 
P  =  0.045  ( s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  *  4,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  29.2  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.83 
Crust2 . 0  =  1.80 
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Station:  GT  BGCA 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0594  (s/km) 
Influence  factor  -  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  94,  G2.50  =  97 
Crustal  Thickness 
HK  method  =  41.9  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.71 
Crust2 . 0  =  1.77 
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Station:  G  ATD 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0547  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  134,  G2.50  =  116 
Crustal  Thickness 
HK  method  =  22.60  km 
Crust2 .0  =  17  km 
Vp/Vs  Ratio 
HK  method  =  1.73 
Crust2 . 0  =  1.79 
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Station:  G  ECH 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0515  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  147,  G2.50  =  130 
Crustal  Thickness 
HK  method  =  25.1  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.72 
Crust2.0  =  1.8 
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Station:  G  GRC 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0464  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  11,  G2.50  =  10 
Crustal  Thickness 
HK  method  =  29.3  km 
Crust2.0  -  31  km 
Vp/Vs  Ratio 
HK  method  -  1.77 
Crust2 . 0  =  1.80 
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Station:  G  MBO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0638  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  12,  G2.50  =  4 
Crustal  Thickness 
HK  method  =  26.3  km 
Crust2.0  =  26  km 
Vp/Vs  Ratio 
HK  method  =  1.98 
Crust2 .0  =  1.84 
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Station:  G  SSB 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0 . 052  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  211,  G2.50  =  192 
Crustal  Thickness 
HK  method  =  271.  km 
Crust2.0  =  31  km 
vp/Vs  Ratio 
HK  method  =  1.69 
Crust2 . 0  =1.8 
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Station:  G  TAM 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0559  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  100,  G2.50  =  30 
Crustal  Thickness 
HK  method  =  40  km 
Crust2.0  -  40  km 
Vp/Vs  Ratio 
HK  method  =  1.6 
Crust2.0  =  1.77 
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Station:  G  WUS 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0648  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  423,  G2.50  =  242 
Crustal  Thickness 
HK  method  =  54.1  km 
Crust2.0  =  50  km 
Vp/Vs  Ratio 
HK  method  =  1.71 
Crust2 . 0  =  1.76 
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Station:  II  AAK 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.065  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  96,  G2.50  =  114 
Crustal  Thickness 
HK  method  =  50  km 
Crust2.0  =  50  km 
Vp/Vs  Ratio 
HK  method  =  1.78 
Crust2.0  =  1.77 
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Station:  II  ABKT 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0573  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  193,  G2.50  =  86 
Crustal  Thickness 
HK  method  =  41.8  km 
Crust2.0  =  46  km 
Vp/Vs  Ratio 
HK  method  =  1.7 
Crust2 . 0  =  1.79 
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Station:  II  BFO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.054  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  -  74,  G2.50  =  62 
Crustal  Thickness 
HK  method  =  26.1  km 
Crust2.0  =  27.5  km 
Vp/Vs  Ratio 
HK  method  =  1.66 
Crust2 . 0  -  1.8 
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Station:  II  ESK 


Shear  Velocity  Models 


Back,  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0625  <s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  11,  G2.50  =  10 
Crustal  Thickness 
HK  method  =  28.2  km 
Crust2 .0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.82 
Crust2 . 0  =  1.79 
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Station:  II  GAR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.065  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  14,  G2.50  =  14 
Crustal  Thickness 
HK  method  =  54.5  km 
Crust2.0  =  50  km 
Vp/Vs  Ratio 
HK  method  =  1.75 
Crust2 . 0  =  1.77 
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Station:  II  KIV 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0551  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  -  0.5 
Number  of  RFTN 
G1.00  =  319,  G2.50  =  408 
Crustal  Thickness 
HK  method  =  1.74  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  = 

Crust2 .0  =  1.78 
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Station:  II  MBAR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0583  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  10,  G2.50  =  9 
Crustal  Thickness 
HK  method  =  34.7  km 
Crust2.0  =  39  km 
Vp/Vs  Ratio 
HK  method  =  1.55 
Crust2 . 0  =  1.77 
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Station:  II  MSEY 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0593  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  10,  G2.50  =  15 
Crustal  Thickness 
HK  method  =  km 
Crust2.0  =  18.7  km 
Vp/Vs  Ratio 
HK  method  = 

Crust2 .0  =  1.79 
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Station:  II  NIL 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0625  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  70,  G2.50  =  42 
Crustal  Thickness 
HK  method  =  51  km 
Crust2.0  =  50  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 .0  =  1.78 
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Group  Velocity  (km/s)  Receiver  Function  (1/s) 
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Station:  II  OBN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0569  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  244,  G2.50  =  206 
Crustal  Thickness 
HK  method  =  50  km 
Crust2.Q  -  46.5  km 
Vp/Vs  Ratio 
HK  method  =  1.7 
Crust2 . 0  =  1.79 
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Station:  II  RAYN 


Shear  Velocity  Models 


Back  azimuth  range  -  All 
Average  ray  parameter 
P  =  0.054  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  72,  G2.50  =  119 
Crustal  Thickness 
HK  method  =  38  km 
Crust2.0  -  41  km 
Vp/Vs  Ratio 
HK  method  =  1.77 
Crust2 . 0  =1.79 
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Station:  IL  IRl 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0-0518  { s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  51,  G2.50  =  52 
Crustal  Thickness 
HK  method  =  48  km 
Crust2 .0=41  km 
Vp/Vs  Ratio 
HK  method  =  1.77 
Crust2 . 0  =  1.786 
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Station:  IL  IR2 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0524  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  31,  G2.50  =  27 
Crustal  Thickness 
HK  method  =  43.2  km 
Crust2.0  =  4 1  km 
Vp/Vs  Ratio 
HK  method  =  1.82 
Crust2 . 0  =  1.786 
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Station:  IL  IR3 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0518  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  56,  G2.50  =  52 
Crustal  Thickness 
HK  method  =  45.2  km 
Crust2.0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.8 
Crust2 .0  =  1.786 
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Station:  IL  IR4 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0509  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  14,  G2.50  =  16 
Crustal  Thickness 
HK  method  =  42  km 
Crust2 .0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.96 
Crust2 . 0  =  1.786 
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Station:  IL  IR5 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0529  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  9,  G2.50  =  12 
Crustal  Thickness 
HK  method  =  49  km 
Crust2 . 0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.85 
Crust2 . 0  =  1.786 
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Station:  IL  IR6 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  -  0.0504  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  7,  G2.50 
Crustal  Thickness 
HK  method  =  44.7  km 
Crust2 . 0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.87 
Crust2 .0  =  1.786 
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Station:  IL_IR7 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0515  (s/km) 
Influence  factor  -  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  32,  G2.50  =  35 
Crustal  Thickness 
HK  method  =  41.5  km 
Crust2.0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.88 
Crust2 .0  =  1.786 


Shear  Velocity  Models 
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Station:  IU  ANTO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.052  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  50,  G2.50  =  33 
Crustal  Thickness 
HK  method  =  35  km 
Crust2.0  =  46  km 
Vp/Vs  Ratio 
HK  method  =  1.71 
Crust2 .0  =  1.77 
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Station:  IU  FURI 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0589  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  12,  G2.50  =  8 
Crustal  Thickness 
HK  method  =  38.60  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.86 
Crust2 .0  =  1.77 
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Station:  IU  GNI 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0571  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  37,  G2.50  =  8 
Crustal  Thickness 
HK  method  =  44.1  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.71 
Crust2 .0  =  1.78 
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Group  Velocity  (km/s)  Receiver  Function  (1/s) 
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Station:  IU  GRFO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0547  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  29,  G2.50  =  32 
Crustal  Thickness 
HK  method  =  29.7  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 . 0  =  1.8 
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Station:  IU  KIEV 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0563  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  81,  G2.50  =  145 
Crustal  Thickness 
HK  method  =  39.5  km 
Crust2.0  =  4 1  km 
Vp/Vs  Ratio 
HK  method  =1.79 
Crust2.0  =  1.78 
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Station:  IU  KMBO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0572  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  1 8 ,  G2.50  =  32 
Crustal  Thickness 
HK  method  -  39.9  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 . 0  =  1.77 
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Station:  IU  KONO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.053  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  13,  G2.50  =  9 
Crustal  Thickness 
HK  method  =  32.2  km 
Crust2 . 0  =  32  km 
Vp/Vs  Ratio 
HK  method  =  1.71 
Crust2 .0  =  1.82 
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Station:  IU  KOWA 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0525  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  2,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  42.6  km 
Crust 2.0  =  km 
Vp/Vs  Ratio 
HK  method  =  1.66 
Crust2.0  = 
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Station:  IU  MSKU 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0498  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  3,  G2.50  =  3 
Crustal  Thickness 
HK  method  =  45.6  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.78 
Crust2 . 0  =  1.79 
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Station:  IU  PAB 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.051  (s/km) 
Influence,  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  44,  G2.50  =  49 
Crustal  Thickness 
HK  method  =  31.5  km 
Crust2.G  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.71 
Crust2 . 0  =1.79 
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Station :  IUTOL 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0548  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  6,  G2.50  =  3 
Crustal  Thickness 
HK  method  =  33.2  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.73 
Crust2 . 0  =  1.79 
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Station:  KN  AML 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0602  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  156,  G2.50  =  107 
Crustal  Thickness 
HK  method  =  53.3  km 
Crust2 .0  =  45  km 
Vp/Vs  Ratio 
HK  method  =  1.93 
Crust2 . 0  =  1.77 
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Station:  KZ  AKT 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0579  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  5,  G2.50  =  12 
Crustal  Thickness 
HK  method  =  44.6  km 
Crust.2.0  =  50  km 
vp/Vs  Ratio 
HK  method  -  1.97 
Crust2 .0  -  1.76 
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Station:  KZ  AKTK 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0608  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  22,  G2.50  =  33 
Crustal  Thickness 
HK  method  =  50.5  km 
Crust2.0  =  50  km 
Vp/Vs  Ratio 
HK  method  =  1.86 
Crust2 .0  =  1.76 
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station:  KZBRVK 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0613  (s/km) 

Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  156,  G2.50 
Crustal  Thickness 
HK  method  =  47.2  km 
Crust2.0  =  4 1  km 
Vp/Vs  Ratio 
HK  method  =  1.73 
Crust2 . 0  =  1.78 


Shear  Velocity  Models 
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Station:  KZCHK 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0623  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  29,  G2.50  =  101 
Crustal  Thickness 
HK  method  =  41.2  km 
Crust2 .0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.72 
Crust2.0  =  1.78 


Shear  Velocity  Models 
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Station:  KZ  KKAR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0626  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  57,  G2.50  =  56 
Crustal  Thickness 
HK  method  =  40.2  km 
Crust2 . 0  =  45  km 
Vp/Vs  Ratio 
HK  method  =  1.65 
Crust2 .0  -  1.77 
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Station:  KZ  PDG 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0642  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  14,  G2.50  =  17 
Crustal  Thickness 
HK  method  =  44.4  km 
Crus t2 .0  =  50  km 
Vp/Vs  Ratio 
HK  method  =  1.86 
Crust2 . 0  =  1.77 
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Station:  KZ  TLG 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.066  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  49,  G2.50  =  41 
Crustal  Thickness 
HK  method  =  46.4  km 
Crust2.0  =  50  km 
Vp/Vs  Ratio 
HK  method  =  1.72 
Crust2.0  =  1.77 
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Station:  KZ  VOS 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0619  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  181,  G2.50  =  183 
Crustal  Thickness 
HK  method  =  38.2  km 
Crust2.0  =  4 1  km 
Vp/Vs  Ratio 
HK  method  =  2.1 
Crust2 . 0  =  1.78 
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Station:  KZ  ZRN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0608  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  38,  G2.50  -  37 
Crustal  Thickness 
HK  method  =  44.5  km 
Crust2 . 0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.74 
Crust2.0  *  1.79 
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Station:  MN  AIO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0445  <s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  2,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  31.2  km 
Crust2.0  =  26.5  km 
Vp/Vs  Ratio 
HK  method  =  1.99 
Crust2.0  =  1.86 
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Station:  MN_AQU 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0586  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  90,  G2.50  =  24 
Crustal  Thickness 
HK  method  =  37.5  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.7 
Crust2 . 0  =  1.79 
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Station:  MNBGY 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0587  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  20,  G2.50  =  7.5 
Crustal  Thickness 
HK  method  =  27.1  km 
Crust2.0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.78 
Crust2 . 0  =1.79 
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Station:  MN  BN I 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0523  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  26,  G2.50  =  15 
Crustal  Thickness 
HK  method  =  45  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.86 
Crust2.0  =  1.77 
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Station:  MN  CEL 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0-0435  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  2,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  27  km 
Crust2.0  =  24.5  km 
Vp/Vs  Ratio 
HK  method  =  1.83 
Crust2 .0  =  1.86 


Time  (s) 


VS  (km/s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 


APPENDIX  B.  INVERSION  SUMMARIES 


156 


Station:  MN  CII 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.041  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  7,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  32  km 
Crust2 .0  =  31  km 
Vp/Vs  Ratio 
HK  method  =  1.84 
Crust2.0  =  1.8 
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Station:  MN  CLTB 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0512  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  5,  G2.50  =  6 
Crustal  Thickness 
HK  method  =  30.1  km 
Crust2.0  =  26  km 
Vp/Vs  Ratio 
HK  method  =  1.78 
Crust2 .0  =  1.83 
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Station:  MN  CUC 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0 . 056  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  3,  G2.50  =  4 
Crustal  Thickness 
HK  method  =  17.9  km 
Crust2 . 0  =  24.5  km 
Vp/Vs  Ratio 
HK  method  =  1.91 
Crust2.0  =  1.86 
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Station:  MN_GFA 

Back  azimuth  range  =  All 


Shear  Velocity  Models 


Average  ray  parameter 

P  =  0.0515  (s/km) 
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Station:  MN  IDI 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0 .  0486  ( s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  -  17,  G2.50  =  5 
Crustal  Thickness 
HK  method  =  28.8  km 
Crust2.0  =  27  km 
Vp/Vs  Ratio 
HK  method  =  1.68 
Crust2 . 0  =1.86 
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Station:  MN  KEG 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0561  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  44,  G2.50  =  44 
Crustal  Thickness 
HK  method  =  32.2  km 
Crust2.0  -  32  km 
Vp/Vs  Ratio 
HK  method  =  1.64 
Crust2 . 0  =  1.81 
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Station:  MN  MDT 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0558  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  19,  G2.50  =  16 
Crustal  Thickness 
HK  method  =  33.2  km 
Crust2.0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.72 
Crust2.0  =  1.77 
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Station:  MN  MEB 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0677  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  -  10,  G2.50  =  3 
Crustal  Thickness 
HK  method  =  26.2  km 
Crust2.0  =  24.5  km 
Vp/Vs  Ratio 
HK  method  =  1.9 
Crust2 .0  =  1.86 
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Station:  MN  TRI 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0532  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  55 /  G2.50  =  41 
Crustal  Thickness 
HK  method  =  36  km 
Crust2.0  ~  37  km 
Vp/Vs  Ratio 
HK  method  =  1.77 
Crust2 . 0  =  1.8 
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Station:  MN  TUE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0-043  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  2,  G2.50  =  1 
Crustal  Thickness 
HK  method  =  30.1  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.88 
Crust2 . 0  =  1.77 
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Station:  MN  VAE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.044  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  1,  G2.50  =  1 
Crustal  Thickness 
HK  method  =  h  km 
Crust2.0  =  26.5  km 
Vp/Vs  Ratio 
HK  method  =  d(h) 

Crust2 . 0  =1.86 
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Station:  MN  VSL 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0516  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  125,  G2.50  =  149 
Crustal  Thickness 
HK  method  =  30.7  km 
Crust2 . 0  =  22  km 
Vp/Vs  Ratio 
HK  method  =  1.69 
Crust2 . 0  =1.82 
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Station:  MN  VTS 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0542  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  73,  G2.50  =  109 
Crustal  Thickness 
HK  method  =  38.1  km 
Crust2.0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.73 
Crust2 . 0  =  1.79 
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Station:  MN  WDD 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0533  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  46,  G2.50  =  33.5 
Crustal  Thickness 
HK  method  =  35.7  km 
Crust2.0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.69 
Crust2 . 0  =  1.81 
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Station:  SR  MAIO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0549  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  8  4 ,  G2.50  =  78 
Crustal  Thickness 
HK  method  =  48.4  km 
Crust2 .0  =  46  km 
Vp/Vs  Ratio 
HK  method  =  1.77 
Crust2.0  =  1.786 
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Station:  XD  AMBA 


Shear  Velocity  Models 


Back  azimuth  range  -  All 
Average  ray  parameter 
P  =  0.0463  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  10,  G2.50  =  7 
Crustal  Thickness 
HK  method  =  34.5  km 
Crust2.0  =  33  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 . 0  =  1.77 
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Station:  XD_BASO 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0558  (s/km) 

Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  16,  G2.50  =  11 
Crustal  Thickness 
HK  method  =  38.2  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.67 
Crust2 . 0  -  1.77 
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Station:  XD  GOMA 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0508  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  25,  G2.50  =  18 
Crustal  Thickness 
HK  method  -  43.3  km 
Crust2.0  -  43.5  km 
Vp/Vs  Ratio 
HK  method  =  1.78 
Crust2 . 0  =  1.77 
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Station:  XD  HALE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0495  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0-5 
Number  of  RFTN 
G1.00  =  21,  G2.50  =  20 
Crustal  Thickness 
HK  method  =  35.2  km 
Crust2.0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.7 
Crust2 .0  =  1.80 
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Station:  XD  INZA 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0526  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  11,  G2.50  =  5 
Crustal  Thickness 
HK  method  =  28.9  km 
Crust2 . 0  =  41.5  km 
Vp/Vs  Ratio 
HK  method  =  1.77 
Crust2 . 0  =1.76 
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Station:  XD_KIBA 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.051  (s/km) 

Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  31,  G2.50  =  23 
Crustal  Thickness 
HK  method  =  39.2  km 
Crust2.0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.72 
Crust2 .0  =  1.77 


Shear  Velocity  Models 
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Station:  XD  KIBE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0505  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  24,  G2.50  =  20 
Crustal  Thickness 
HK  method  =  38.1  km 
Crust2.0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.7 
Crust2 . 0  =  1.77 
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Station:  XD  KOMO 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.052  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  30 ,  G2.50  =  22 
Crustal  Thickness 
HK  method  =  37  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.71 
Crust2 .0  =  1.77 
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Station:  XDKOND 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0494  (s/km) 

Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  18,  G2.50  =  17 
Crustal  Thickness 
HK  method  =  38.3  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.66 
Crust2 .0  =  1.77 
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Station:  XDLONG 

Back  azimuth  range  =  All 


Shear  Velocity  Models 


Average  ray  parameter 

P  =  0.0521  (s/km) 
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Station:  XD  MB WE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0513  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  19,  G2.50  =  15 
Crustal  Thickness 
HK  method  =  36.5  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.75 
Crust2 . 0  =  1.77 
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Station:  XD  MITU 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0532  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  20,  G2.50  =  15 
Crustal  Thickness 
HK  method  =  37.2  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  -  1.74 
Crust2.0  =  1.77 
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Station:  XD  MTAN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.049  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  4,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  37.9  km 
Crust2 . 0  =  41.5  km 
Vp/Vs  Ratio 
HK  method  =  1.67 
Crust2 . 0  =  1.76 
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Station:  XD  MTOR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.051  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  38,  G2.50  =  30 
Crustal  Thickness 
HK  method  =  37.9  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.74 
Crust2 . 0  =  1.77 
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Station:  XD  PAND 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0472  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  1 6 ,  G2.50  =  11 
Crustal  Thickness 
HK  method  =  32.4  km 
Crust2 .0  =  33  km 
Vp/Vs  Ratio 
HK  method  =  1.8 
Crust2 .0  =  1.77 
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Station:  XD  PUGE 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0527  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  20,  G2.50  =  16 
Crustal  Thickness 
HK  method  =  35.7  km 
Crust2.0  =  37  km 
Vp/Vs  Ratio 
HK  method  =  1.77 
Crust2 .0  =  1.76 
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Station:  XD  RUNG 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P' =  0.053  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Humber  of  RFTN 
G1.00  =  25 ,  G2.50  =  20 
Crustal  Thickness 
HK  method  =  39.2  km 
Crust2.0  =  41.5  km 
Vp/Vs  Ratio 
HK  method  =  1.79 
Crust2 .0  =  1.76 
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Station:  XD  SING 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0537  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  -  0.5 
Number  of  RFTN 
G1.00  =  24,  G2.50  =  19 
Crustal  Thickness 
HK  method  =  38.5  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.78 
Crust2.0  =  1.77 
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Station:  XD  TARA 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0499  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  25,  G2.50  =  21 
Crustal  Thickness 
HK  method  =  37.2  km 
Crust2 . 0  -  38  km 
Vp/Vs  Ratio 
HK  method  =  1.75 
Crust2 .0  =  1.77 
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Station:  XD  TUND 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0531  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  22,  G2.50  =  15 
Crustal  Thickness 
HK  method  -  34.6  km 
Crust2.0  =  33  km 
Vp/Vs  Ratio 
HK  method  =  1.7 
Crust2 . 0  =  1.77 
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Station:  XD  URAM 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0556  (s/km) 
Influence  factor  -  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  21,  G2.50  =  15 
Crustal  Thickness 
HK  method  =  40.3  km 
Crust2 .0  =  37  km 
Vp/Vs  Ratio 
HK  method  =  2 
Crust2.0  =  1.76 
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Station:  XG  AGIN 


Shear  Velocity  Models 


Back  azimuth  range  =  Ail 
Average  ray  parameter 
P  =  0.0527  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  21,  G2.50  =  4 
Crustal  Thickness 
HK  method  =  37.3  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 .0  =  1.77 
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Station:  XG  AHLT 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0539  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  22,  G2.50  =  16 
Crustal  Thickness 
HK  method  =  42.9  km 
Crust2 . 0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 .0  =  1.77 
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Station:  XG  BNGL 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0521  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  29,  G2.50  =  32 
Crustal  Thickness 
HK  method  =  42.8  km 
Crust2 . 0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.84 
Crust2 . 0  =  1.77 
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Station:  XG  BTLS 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0528  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  26,  G2.50  =  27 
Crustal  Thickness 
HK  method  =  40.5  km 
Crust2 .0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.83 
Crust2 . 0  =1.77 


-5  0  5  10  15  20  25  30  35  40  45  50  55  60 

Time  (s) 


VS  (km/s) 


Period  (s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 


APPENDIX  B.  INVERSION  SUMMARIES 


196 


Station:  XG  BYBT 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0514  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  20,  G2.50  =  13 
Crustal  Thickness 
HK  method  =  38.8  km 
Crust2.0  -  36  km 
Vp/Vs  Ratio 
HK  method  =  1.89 
Crust2 . 0  =  1.79 
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Station:  XG  BYKN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0516  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  31,  G2.50  =  14 
Crustal  Thickness 
HK  method  =  46.2  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.63 
Crust2 . 0  =1.77 
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Station:  XG  CMCY 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0548  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  9 ,  G2.50  =  4 
Crustal  Thickness 
HK  method  =  41.7  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.83 
Crust2 . 0  =  1.77 
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Station:  XG  DGRL 


Shear  Velocity  Models 


Back  azimuth  range  =  Ail 
Average  ray  parameter 
P  =  0.0551  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  19,  G2.50  =  10 
Crustal  Thickness 
HK  method  =  37  km 
Crust2.0  =  39  km 
Vp/Vs  Ratio 
HK  method  =  1.86 
Crust2 . 0  =  1.76 
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Station:  XG  DGSU 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0535  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  13,  G2.50  =  6 
Crustal  Thickness 
HK  method  =  4 1  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.84 
Crust2 . 0  =  1.77 


-5  0  5  10  15  20  25  30  35  40  45  50  55  60 


Time  (s) 


VS  (km/s) 


Period  (s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 


APPENDIX  B.  INVERSION  SUMMARIES 


201 


Station:  XG  DYBR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0539  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G 1 . 00  =  22,  G2.50  =  18 
Crustal  Thickness 
HK  method  =  35.3  km 
Crust2 . 0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.67 
Crust2.0  =  1.78 
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Station:  XG  ERGN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0529  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  22,  G2.50  =  18 
Crustal  Thickness 
HK  method  =  45.2  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.7 
Crust2 . 0  =  1.77 
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Station:  XG  EZRM 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0542  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  19,  G2.50  =  4 
Crustal  Thickness 
HK  method  =  40.3  km 
Crust2 .0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.79 
Crust2 . 0  =  1.79 
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Station:  XG  HAMR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0-047  ( s/kxn) 
Influence  factor  =0-5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  12,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  45.6  km 
Crust2 . 0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.69 
Crust2 . 0  =  1.77 
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Station:  XG  HINS 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0504  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  21,  G2.50  =  25 
Crustal  Thickness 
HK  method  =  44.2  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.57 
Crust2 . 0  =  1.77 


Time  (s) 


VS  (km/s) 


Period  (s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 


APPENDIX  B.  INVERSION  SUMMARIES 


206 


Station:  XG  HRPT 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0491  (s/km) 
Influence  factor  -  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  20,  G2.50  =  38 
Crustal  Thickness 
HK  method  =  35.6  km 
Crust2.Q  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.81 
Crust2 .0  =  1.77 
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Station:  XG  HRSN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0515  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  14,  G2.50  =  6 
Crustal  Thickness 
HK  method  =  41  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.81 
Crust2.0  =  1.77 
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Station:  XG  HUB03 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0529  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  19,  G2.50  =  10 
Crustal  Thickness 
HK  method  =  50  km 
Crust2.0  =  65  km 
Vp/Vs  Ratio 
HK  method  =  1.66 
Crust2.0  =  1.76 
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Station:  XG  ILIC 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0529  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  40,  G2.50  =  24 
Crustal  Thickness 
HK  method  =  42  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.91 
Crust2 . 0  =1.77 
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Station:  XG  IMRL 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.057  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  5,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  43.3  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.84 
Crust2 . 0  =1.77 
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Station:  XG  KARS 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0498  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  12,  G2.50  =  5 
Crustal  Thickness 
HK  method  =  36.5  km 
Crust2.0  =  39  km 
Vp/Vs  Ratio 
HK  method  =  2.17 
Crust2 . 0  =  1.76 
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Station:  XG  KOTK 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0527  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  28,  G2.50  =  13 
Crustal  Thickness 
HK  method  =  40.85  km 
Crust2.0  =  39  km 
Vp/Vs  Ratio 
HK  method  =  1.83 
Crust2.0  =  1.76 


0.8 


0.6  - 


0.4 


0.2  - 


0.0  - 


Observed 

Harvard 

Pasyanos 


1  I  "  "  !  1  "  I  "  "  I  1 


-5  0  5  10  15  20  25  30  35  40  45  50  55  60 

Time  (s) 


VS  (km/s) 


Period  (s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 


APPENDIX  B.  INVERSION  SUMMARIES 


213 


Station:  XG  KRLV 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0516  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  24  ,  G2.50  =  11 
Crustal  Thickness 
HK  method  =  40.4  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.83 
Crust2 .0  =  1.77 
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Station:  XG  KTLN 


Shear  Velocity  Models 


Back  azimuth  range  -  All 
Average  ray  parameter 
P  =  0.0509  (s/km) 
Influence  factor  -  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  10,  G2.50  =  11 
Crustal  Thickness 
HK  method  =  24.7  km 
Crust2.0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  2.19 
Crust2 . 0  =  1.78 
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station:  XG  KYPR 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0515  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  36,  G2.50  =  31 
Crustal  Thickness 
HK  method  =  32.6  km 
Crust2.0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.85 
Crust2 .0  =  1.78 
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Station:  XG  MRDN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0512  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  25 ,  G2.50  =  19 
Crustal  Thickness 
HK  method  =  35.9  km 
Crust2.0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2.0  =  1.78 
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Station:  XG  MSDY 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  -  0.0511  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  19,  G2.50  =  23 
Crustal  Thickness 
HK  method  =  38.1  km 
Crust2 . 0  =  4 1  km 
vp/vs  Ratio 
HK  method  =  1.84 
Crust2 . 0  =  1.79 
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Station:  XG_MUSH 

Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0505  {s/km) 

Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  21,  G2.50  =  17 
Crustal  Thickness 
HK  method  =  44.7  km 
Crust2-0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.75 
Crust2 . 0  =  1.77 


Shear  Velocity  Models 
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Station:  XG  SILN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0511  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  12,  G2.50  =  12 
Crustal  Thickness 
HK  method  =  41.3  km 
Crust2.0  =  43  km 

Vp/Vs  Ratio 

HK  method  =  1.64 
Crust2 . 0  =  1.77 
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Station:  XG  SIRN 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0535  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  30,  G2.50  =  13 
Crustal  Thickness 
HK  method  =  38.5  km 
Crust2.0  =  38  km 
Vp/Vs  Ratio 
HK  method  =  1.68 
Crust2 . 0  =1.8 
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Station:  XG  TRB01 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.063  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  9,  G2.50  =  12 
Crustal  Thickness 
HK  method  =  71.6  km 
Crust2.0  =  65  km 
vp/vs  Ratio 
HK  method  =  1.69 
Crust2.0  =1.76 


Time  (s) 


VS  (km/s) 


Group  Velocity  (km/s)  Receiver  Function  (1/s) 


APPENDIX  B.  INVERSION  SUMMARIES 


222 


Station:  XG  TTB01 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0665  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  2,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  68.84  km 
Crust2.0  =  65  km 
Vp/Vs  Ratio 
HK  method  =  1.68 
Crust2 . 0  =  1.76 
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Station:  XG  UZML 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  -  0.049  { s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  12,  G2.50  =  28 
Crustal  Thickness 
HK  method  =  42.4  km 
Crust2.0  =  43  km 

Vp/Vs  Ratio 

HK  method  =  1.79 
Crust2 .0  =  1.77 
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Stations  XI  AFIF 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.049  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  23,  G2.50  =  29 
Crustal  Thickness 
HK  method  =  36.35  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.75 
Crust2 .0  =  1.77 
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Station:  XI  BISH 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.043  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  2,  G2.50  =  2 
Crustal  Thickness 
HK  method  =  37.13  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.82 
Crust2 . 0  =  1.77 
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Station:  XI  HALM 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0485  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  31,  G2.50  =  31 
Crustal  Thickness 
HK  method  =  34.6  km 
Crust2.0  -  43  km 
Vp/Vs  Ratio 
HK  method  =  1.79 
Crust2 . 0  =1.79 
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Station:  XI  RANI 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0501  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  18,  G2.50  =  18 
Crustal  Thickness 
HK  method  =  38.1  km 
Crust2 . 0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 . 0  =  1.77 
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Station:  XI  RIYD 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0557  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  13,  G2.50  =  9 
Crustal  Thickness 
HK  method  =  41.6  km 
Crust2.0  =  41  km 
Vp/Vs  Ratio 
HK  method  =  1.64 
Crust2 . 0  =  1.79 
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Station:  XI  SODA 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.053  (s/km) 
Influence  factor  =0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  31,  G2.50  =  33 
Crustal  Thickness 
HK  method  =  38.4  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.73 
Crust2 . 0  =  1.77 
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Group  Velocity  (km/s)  Receiver  Function  (1/s) 
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Station:  XI  TAIF 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.045  ( s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1 .00  =  4,  G2.50  =  4 
Crustal  Thickness 
HK  method  =  37.4  km 
Crust2.0  =  36  km 
Vp/Vs  Ratio 
HK  method  =  1.74 
Crust2 . 0  =  1.77 
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Station:  XI_UQSK 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0496  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  11,  G2.50  =  11 
Crustal  Thickness 
HK  method  =  36.5  km 
Crust2.0  =  43  km 
Vp/Vs  Ratio 
HK  method  =  1.74 
Crust2 . 0  =  1.77 
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Station:  XWAHQI 


Shear  Velocity  Models 


Back  azimuth  range  =  All 
Average  ray  parameter 
P  =  0.0641  (s/km) 
Influence  factor  =  0.5 
Smoothing  factor  =  0.5 
Number  of  RFTN 
G1.00  =  29,  G2.50  =  13 
Crustal  Thickness 
HK  method  =  57.3  km 
Crust2.0  =  50  km 
Vp/Vs  Ratio 
HK  method  =  1.76 
Crust2 .0  =  1.76 
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